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Abstract 

\ 

This  report  studies  three  tolman  Filter  implementations  for  the  de- 
orbit/reentry phase  of  a Space  Shuttle  Crbiter  mission  profile.  Each 
design  uses  measurements  from  Global  Positioning  System  satellites  to 
update  Kalman  Filter  estimates  of  selected  error  states  and  to  bound  the 
Shuttle's  INS  error  grcvrth.  Covariance  analysis  techniques  are  used  to 
compare  the  performances  of  17-state,  11-state,  and  S-state  filter  de- 
signs using  a digital  computer  simulation  of  a NASA  CFT-1  Space  Shuttle 
mission  trajectory.  A 73-state  truth  model  of  the  Space  Shuttle  II-IU, 
the  Crbiter  GPS  receiver  Master  Tuning  Unit,  and  the  GPS  satellite  con- 
stellation is  used  as  a system  simulation  to  examine  the  performance  of 
the  design  filters.  The  results  are  shown  primarily  in  the  form  of 
computer-generated  plots  of  both  filter  indicated  errors  and  true  system 
errors  versus  time.  It  is  concluded  from  the  results  that  an  8-state 
filter  design  is  able  to  perform  well  enough  to  be  a primary  candidate 
for  design  implementation,  and  also  that  further  refinements  are  neces- 
sary in  the  user  clock  model  of  the  design  Kalman  Filter.  ^ 


A K/\LMAN  FILTER  LE  -IGN  FOR  THE  SPACE  SHUTTLE 
CRBITER  INERTIAL  MEASURING  UNIT  DURIN’G  LECRBIT/REENTRY 
USING  GLOBAL  POSITIONING  SYSTEM  SATELLITE  INFORI-IATION 


I.  Introduction 


Purpose 

This  study  is  an  investigation  of  mininum  Kalman  Filter  state  size 
required  to  achieve  specified  accuracy  from  the  Space  Shuttle  Orbiter 
(hereafter  referred  to  as  the  Orbiter)  aided  inertial  navigation  system 
(IKS)  during  deorbit/reentry  operation.  The  Kalman  Filter  designs  stud- 
ied utilize  range  and  pseudo  range-rate  measurements  from  the  Air  Force's 
Glo’oal  Positioning  System  (GPS)  satellites  as  an  updating  aid  for  bound- 
ing the  INS's  inherent  error  growth, 

NASA's  Orbiter  uses  three  KT-70  based  Inertial  Measuring  Units 
(iMUs)  for  its  navigation  and  for  providing  much  of  its  guidance  infor- 
mation during  de orbit/reentry,  NASA  proposes  to  use  a Kalman  Filter,  a 
recursive  data  processing  algorithm  residing  in  the  Orbiter 's  onboard 
computer,  to  limit  li'IU  error  growth  during  the  deorbit/reentry  portion 
of  the  mission.  The  proposed  filter  (Ref  10j2.7  and  uses  a 

pseudo  drag  measurement  scheme  to  compute  a drag  update  measurement  for 
the  Kalman  Filter  during  initial  deceleration  and  ionospheric  blackout 
(Ref  9*^*^  3-nd  IOjA.I).  This  involves  computation  of  a pseudo  measure- 
ment of  drag  from  the  IMU  accelerometer  outputs  and  comparison  of  these 
measurements  to  expected  values  of  drag  generated  from  current  estimates 
of  position  and  velocity  from  the  LMU  and  from  table  look  up  routines. 
This  was  to  be  the  sole  updating  aid  for  the  IMU  until  after  exit  from 
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ionospheric  radio  blackout.  At  this  time  TAGAN  and  baro  altimeter  meas- 
urements would  replace  pseudo  drag  as  the  updating  aids  for  the  Kalman 
Filter. 

This  study  investigates  the  use  of  GPS  measurements  as  updating 
aids  for  the  filter  from  deorbit  ignition  onward  to  ionospheric  black- 
out f and  again  after  blackout  ends  rather  than  TAGAN  and  baro  altimeter 
measurements.  Measurements  would  be  unavailable  for  a brief  period 
during  ionospheric  blackout,  during  which  time  pseudo  drag  would  be  used 
to  bound  the  rate  of  error  growth.  After  blackout  ends  GPS  measurements 
would  continue  until  the  Orbiter  reached  20,000  feet  above  mean  sea. 
level  (MSL)  v/hen  the  microwave  landing  system  (MLS)  ’would  take  over. 
Except  for  the  period  of  ionospheric  blackout,  GPS  measurements  are  the 
sole  updating  aids  used  in  the  filter  designs  studied. 

Deorbit/Reentry  of  the  Space  Shuttle  Orbiter 

The  data  for  the  simulated  trajectory  in  this  study  is  from  a NASA- 
supplied  data  tape  (Ref  16)  for  NASA's  early  Orbiter  mission  OFT-1 
trajectory  (Ref  8j6).  The  deorbit  operation  begins  with  deorbit  igni- 
tion, which  occurs  approximately  over  New  Zealand  for  this  trajectory, 
and  ends  with  landing  at  Edwards  Air  Force  Base,  Galifornia  (see  Fig  l). 
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Deorbit  ignition  occurs  at  an  altitude  of  72? >550  feet.  Ionospheric 


blackout  is  not  well  defined  for  the  L 'oand  GPS  signals,  but  for  this 
study  is  assumed  to  occur  at  an  altitude  of  approximately  293>000  feet, 

22  minutes,  11  seconds  after  deorbit  ignition.  It  is  assumed  to  end  at 
177>500  feet,  10  minutes,  45  seconds  later,  at  v.’hich  time  GPS  signals 
are  again  available  for  processing.  From  this  point  until  the  KLS  takes 
over  there  is  another  9 minutes,  37  seconds  during  which  GPS  measurements 
are  used. 


Global  Positioning  System 

The  Air  Force's  Global  Positioning  System  (GPS)  vjill  eventually 
provide  world vfide  navigation  capability  to  users  that  are  able  to  receive 
signals  transmitted  by  navigation  satellites  currently  being  tested. 

Very  accurate  position  and  velocity  information  will  bo  available  for 
those  users  that  can  process  all  the  information  on  the  GPS  satellite 
signals  and  that  can  select  an  optimum  configuration  of  four  satellites 
out  of  those  that  are  visible  to  them  (Ref  1:18-19). 

The  Space  Shuttle  will  be  able  to  use  this  information  to  maintain 
accurate  position  and  velocity  fixes  without  having  to  resort  to  ground- 
link  information.  This  information  can  also  be  used  to  place  a bound  on 
the  error  growth  inherent  in  the  onboard  H-iU  through  the  use  of  Kalman 
Filter  error  estimation. 

This  study  simulates  a 24  satellite  constellation  tliat  propagates 
in  time  as  actual  satellites  vrould  (see  Chapter  III  and  Appendix  B). 

The  simulation  employs  a satellite  selection  criterion  to  obtain  a 
"best"  four-satellite  configuration  of  the  satellites  that  are  visible 
to  the  Crbiter.  Measurements  from  these  four  satellites  in  the 
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simulation  are  used  to  update  a simulation  of  the  Kalman  Filter  estimate 
of  errors  being  committed  by  the  Crbiter  IKU  and  the  filter  then  applies 
a correction  to  the  IKU  itself  based  on  its  updated  estimate  of  these 
errors. 

Assumptions  and  Limitations 

A fundamental  assumption  of  this  study  is  that  a radio  blackout 
period  of  specified  length  exists  during  reentry  due  to  ionospheric 
effects.  During  this  period  of  blackout,  GPS  signals  are  not  available 
for  processing.  It  is  assumed  that  drag  update,  described  earlier,  will 
be  used  during  this  bla.ckout  period  to  limit  error  grovjth.  Due  to 
limitations  of  time  and  resources,  drag  updating  was  not  implemented 
in  the  simulations  of  this  study.  Because  of  this,  the  performance  of 
the  design  filters  were  investigated  in  two  phases,  from  deorbit  igni- 
tion to  blackout,  and  from  exit  blackout  to  20,000  feet. 

An  initialized  covariance  of  errors  matrix  at  deorbit  ignition  was 
supplied  by  NASA  (Ref  5)  for  position  and  velocity  errors.  Initial 
conditions  on  other  errors  were  taken  from  studies  done  for  NASA  by  The 
Analytic  Sciences  Corporation  (TASG)  (Ref  10:2,38-2,40  and  11:2,13- 
2,14)  as  vfere  some  of  the  initial  conditions  for  the  filter.  These 
initial  conditions  assume  an  uploaded  navigation  fix  and  an  IMU  align- 
ment by  star  tracker  two  hours  prior  to  deorbit,  with  no  GPS  measure- 
ments taken  until  time  for  deorbit,  TASC  supplied  an  initialized 
covariance  matrix  for  the  position  and  velocity  errors  at  exit  black- 
out (Ref  12),  The  values  in  the  matrix  assume  GPS  measurements  until 
blackout  and  the  use  of  drag  updating  during  blackout  (Ref  11:2,23- 
2,24), 
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Even  though  the  Global  Positioning  System  will  not  be  completed 
until  193^>  a completed  system  of  24  satellites  is  assumed.  This  is  be- 
cause the  use  of  test  phase  GPS  configurations  has  been  investigated  by 
TASG  already  (Ref  11: Chap  2)  and  the  aim  of  this  study  was  to  investigate 
Kalman  Filter  design  to  use  GPS  signals  rather  than  to  investigate  the 
effects  of  degraded  GPS  configuration  upon  filter  performance.  Also 
assumed  is  the  use  of  a high  accuracy  receiver  on  the  Orbiter  which  is 
capable  of  processing  all  information  available  on  the  GPS  signals  and 
designed  for  the  type  of  highly  dynamic  platform  that  the  Orbiter  will 
be  (Ref  10:1.6). 

Finally,  it  is  assumed  that  linear  approximations  to  the  real  world 
can  be  made  so  that  linear  differential  equation  description  can  be  used 
to  form  a model  of  the  IMU  error  description.  This  is  not  very  re- 
strictive and  should  not  affect  the  validity  of  the  study.  This  will 
allow  the  Covariance  Analysis  to  be  used  as  a tool  in  designing  the 
l^lman  Filter.  In  forming  the  models  that  describe  the  errors,  Gaussian 
descriptions  and  whiteness  of  noises  used  in  the  models  are  both  assumed 
so  that  the  errors  can  then  be  described  as  Gaussian,  characterized  by 
their  means  and  variances. 

Report  Organization 

This  report  will  present  the  elements  of  the  study  in  the  following 
order  by  chapter.  Chapter  II  will  discuss  Space  Shuttle  coordinate  sys- 
tems and  transformations.  Range  and  range  rate  measurei..  nts  using  GPS 
will  be  presented  in  Chapter  III.  Covariance  propagation  equations  used 
in  the  simulation  are  described  in  Chapter  IV.  In  Chapter  V the  "Truth 
Model"  description  of  the  Space  Shuttle  IMU  and  GPS  is  given,  and 
Chapter  VI  describes  the  Kalman  Filter  designs  that  were  investigated 
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in  this  study.  Chapter  VII  presents  the  results  of  the  study  and  the 
conclusions  drawn  from  these  results.  Appendix  A is  composed  of  the 
computer  generated  plots  that  present  detailed  results  of  the  filter- 
indicated  filter  performance  and  the  true  filter  performance.  Appendix 
B presents  a general  description  of  the  computer  program  used  to  per- 
form the  covariance  analysis  and  a depiction  of  those  portions  of  the 
program  unique  to  this  study.  Appendix  C presents  filter  design  param- 
eters recommended  in  other  studies  and  compares  them  to  the  final  design 
values  of  this  study.  Appendix  D consists  of  an  addendum  presentiTig 
results  for  the  8-state  filter  design  with  a corrected  schedule  of 
specific  force  to  the  truth  model  clock  states.  It  will  present  per- 
formance plots  and  a discussion  of  the  results. 

0 


II,  Space  Shuttle  Coordinate  Systems  and  Transformations 


Reference  Coordinate  Frames 

The  coordinate  frames  of  interest  in  describing  the  motion  of  the 
Orbiter  and  the  measurements  of  the  GPS  range  and  range-rate  are  des- 
cribed in  the  following  figures  and  paragraphs.  These  are  the  coordinate 
frames  specified  by  NASA  (Ref  4)  and  described  and  used  by  TASC  in  NASA 
sponsored  studies  relating  to  Space  Shuttle  deorbit/reentry  navigation 
(Ref  9sA,0-A,5)  from  vdiich  the  following  definitions  come. 


Earth  Centered  Inertial  Frame 

The  coordinate  frame  most  important  to  this  study  was  the  earth 
centered  inertial  frame  or  ECIF,  vjhich  will  be  denoted  by  a superscript 
i.  Of  course,  this  is  only  an  approximation  to  an  inertial  reference 
frame,  being  attached  to  the  earth,  but  since  it  is  nonrotating  with 
respect  to  the  "fixed"  stars,  for  short  term  use  it  is  an  adequate 
approximation.  It  is  in  this  reference  frame  that  the  navigation  calcu- 
lations are  performed,  the  GPS  simulation  is  made,  and  most  of  the 
trajectory  data  is  provided.  It  can  be  described  verbally  as  having 
the  axis  positive  through  the  earth's  spin  axis,  the  axis  through 
Aries  of  1950 » the  axis  completing  the  third  axis  of  the  right 
hand  triad  (see  Figure  2), 

Each  of  the  other  reference  frames  are  related  to  the  ECIF  frame 
by  transformation  matrices,  described  in  the  next  sections.  The  trans- 
formation matrices  either  transform  from  a frame,  x,  to  the  ECIF, 

or  to  a frame,  x,  -from  the  ECIF. 
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East-Korth-Un  Reference  Frame 


The  local  level,  east-north-up  (ENU)  frame  is  used  in  this  study  in 
the  algorithm  that  calculates  vfhich  four  satellites  to  use  out  of  the 
set  of  all  visible  satellites.  This  frame  is  best  described  by  the  way 
in  which  it  is  generated.  The  up  unit  vector,  1^  is  described  as  being 
along  the  Crbiter's  position  vector,  from  earth  center  to  the  vehicle. 
The  east  unit  vector,  IE  is  described  as  the  resultant  of  the  cross 
product  of  the  polar  unit  vector,  12^  and  the  Crbiter  position  vector, 

RV.  The  north  unit  vector  is  then  the  cross  product  of  IJJ  and  IE. 

This  calculation  is  described  mathematically  as  follows i 


RU 


(1) 


IZ^  X m 
X m 


(2) 


IN  = lU  X IE 


(3) 


The  direction  cosine  matrix  to  transform  from  the  EC IF  to  the  ENU  is 
then« 


IE 


(4) 


O O •*  • • • V*  ' * * ' - ' •’f  - V*  • - • v_'  • • • * • • . » i.  . • r * V V r •*  . . 


Relative  Velocity  Fnrr.e 

The  frame  which  is  used  to  relate  gravity  deflections  and  anomaly  to 
the  calculations  is  the  Relative  Velocity  Frame  (v),  related  to  the 
Orbiter's  earth  relative  velocity  vector.  The  unit  vectors  which  are  used 
to  define  the  transform  from  the  V frame  to  the  EGIF  frame  for  computation 
purposes  are  defined  by  the  use  of  the  unit  vector  lU  and  the  position 
vector  defined  previously,  IVR,  the  unit  vector  along  the  Orbiter's 
earth  relative  velocity  vector,  and  W,  the  Orbiter's  inertially  refer- 
enced velocity.  Mathematically,  this  transform  is  described  asi 


IVR  = 


VU  - WIE  X RU 


VU  - WIE  X RU 


where  WIE  is  the  earth's  spin  vector.  Furthermore,  letting* 


IW  ■=  lU  X IVR 


the  transform  is  defined  as* 


^i/V 


t t t 

lU  IVR  IW 

T i T 


Platform  Frame 

The  final  frame  of  interest  is  the  platform  coordinate  frame,  or  P 
frame,  in  which  the  accelerometer  measurements  and  gyro  bias  errors  are 


mo? 
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defined.  The  Crbiter  uses  a space  stable  IMU  platform,  which  is  defined 
at  launch,  and  thus  the  transform  from  the  platform  frame  to  the  EGIF 
frame  is  a simple  rotation  which  remains  constant  throughout  the  mission 
The  transform  direction  cosine  to  go  from  the  P frame  to  the  ECIF  frame 
was  provided  by  NASA  (Ref  16)  and  is  given  by 

5.571156-4-  E-1  -1.4125705  E-1  -8.1333281  E-1 

T^^p  = -1.4473432  E-1  -9.8686101  E-1  7.I8I3144  E-2  (7) 

-8.1772438  E-1  7.8433026  E-2  -5.7-24052  E-1 


III.  Global  Positioning  System 


t! 


Description 

Currently  the  Air  Force  and  Navy  are  jointly  developing  a system 
that  will  use  signals  from  a network  of  satellites  to  provide  worldwide 
navigation  capability  with  high  accuracy  to  suitably  equipped  users. 

At  present  the  system  is  entering  a test  and  evaluation  period  (Phase  l) 
during  which  a limited  number  of  satellites  will  be  used  to  prove  con- 
cepts and  evaluate  initial  designs.  When  complete  (Phase  III)  the  sys- 
tem will  consist  of  three  major  segments  which  are  to  be  called  the 
space  segment,  the  user  segment,  and  the  control  segment.  An  inter- 
mediate period  (Phase  II)  for  further  refinement  of  equipment  and 
evaluation  may  or  may  not  be  implemented  depending  on  need  and  cost. 
Completion  of  the  system  is  presently  scheduled  for  1984  (Ref  3«4.7). 

In  Phase  III  the  space  segment  will  consist  of  24  satellites 
divided  into  groups  of  8 satellites  spread  out  in  each  of  three  evenly 
spaced,  circular,  nonsynchronous  orbits.  The  satellites  will  be  in 
an  orbit  of  period  equal  to  11  hours  58  minutes  so  that  each  satellite 
will  repeat  a ground  track  once  every  24  hours.  The  satellites  will 
use  cesium  standard  clocks  to  provide  a very  accurate  clock  signal 
with  low  drift  characteristics.  Each  satellite  will  transmit  coded 
signals  with  highly  directional  chiaracteristics.  These  signals  will 


y 

y 

<.• 

K." 

k 
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allow  a user  to  determine  his  own  relative  position,  velocity  and 
clock  error  factors  by  using  the  information  coded  into  the  signals. 
The  user  segment  will  consist  of  an  L band  antenna , an  onboard 
clock  or  timing  unit,  a receiver  capable  of  processing  all  or  part 
of  the  signal  available  from  the  satellites,  and  an  onboard  computer 
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to  process  the  information.  When  a user's  set  is  to  be  employed,  it 
will  scan  above  its  local  horizon  and  automatically  select  four  satel- 
lites from  those  visible  that  will  provide  the  best  geometrical  config- 
uration to  give  the  best  accuracy  (see  Appendix  3).  The  measurements 
received  from  each  of  the  four  satellites  selected  will  then  provide  a 
basis  for  the  computer  solution  to  three  components  of  position  plus 
the  user  clock  bias  factor,  and  if  desired,  the  solution  to  three 
components  of  velocity  plus  the  user  clock  frequency  error. 

The  control  segment  will  consist  of  four  monitor  stations  located 
in  the  United  States  and  Guam  and  an  upload  station  at  Vandenberg  Air 
Force  Base  in  California.  These  stations  will  monitor  the  accuracy  of 
each  satellite's  signals  and  the  actual  position  of  each  satellite. 

Once  each  day,  as  a satellite  makes  its  closest  approach  to  Vandenberg, 
updated  clock  and  ephemeris  information  will  be  uploaded  to  the 
satellite. 


Using  GPS 

The  way  in  which  the  satellites  provide  the  information  to  the 
user  will  be  by  transmitting  tvfo  distinct  signals  on  separate  L band 
frequencies,  giving  the  time  the  signal  was  transmitted  plus  the 
satellite  position  and  velocity  at  that  time.  3y  measuring  the 
Doppler  shift  in  the  signals  the  user  will  be  able  to  compute  a range 
rate  to  the  satellite.  By  synchronizing  its  clock  to  those  of  the 
satellites,  the  time  difference  from  signal  transmission  to  reception 
can  be  used  to  determine  range  to  the  satellite.  Use  of  the  two 
separate  L band  signals  will  allow  the  user  to  determine  the  varia- 
tion in  velocity  computed  due  to  ionospheric  effects  on  the  signals. 


The  fact  that  the  satellite  will  be  uploaded  daily  is  an  indication 


that  there  is  some  decay  in  the  accuracy  of  information  to  be  expected 
as  well  as  some  slight  clock  drift.  The  time  since  uploading  will  thus 
be  a factor  in  the  accuracy  of  the  information  obtained  by  the  user, 
and  so  also  in  the  magnitude  of  the  error  committed  by  using  this  infor- 
mation. Typical  magnitudes  of  RKiS  range  error  are  approximately  five 
feet  just  prior  to  upload  and  approximately  one  foot  just  after  upload. 
Comparable  magnitudes  in  the  RMS  range  rate  errors  are  approximately 
0.003  feet  per  second  just  prior  to  upload  and  O.OOO5  feet  per  second 
just  after  upload  (Ref  10*2.8). 
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IV,  Error  Covariance  Propagation  Equations 


Covariance  Analysis  Equations 

This  chapter  presents  the  basic  covariance  time  propagation  equa- 
tions, The  equations  propagate  the  filter  estimated  error  covariance 
and  the  system  actual  error  covariance.  The  notation  used  will  also  be 
presented.  The  derivation  of  the  Kalman  Filter  is  well  documented  in 
the  literature  and  the  references  provide  two  sources  (Ref  13* 5* 7“5*^0 
and  14« 247-261),  The  actual  implementation  within  the  General 
Covariance  Analysis  Program  (GCAP)  (Ref  6)  is  presented  in  Appendix  3, 
The  error  model  equations  for  the  system,  or  "truth"  model,  and  for  the 
suboptimal  Kalman  Filter  are  presented  in  the  next  chapters. 

The  Kalman  Filter  requires  a linear  or  linearized  differential 
equation  description  of  the  system.  This  is  acknowledged  to  be  only  an 
approximation  to  the  real  world  nonlinear  model  which  would  be  much 
more  difficult  to  implement  on  a digital  computer.  The  system  model 
used  is  an  integrated  simulation  of  all  of  the  significant  errors 
committed  by  the  Space  Shuttle  H'iU  and  by  GPS,  These  srroi's  include 
both  those  errors  caused  by  imperfections  in  the  system  sensors,  such 
as  gyro  drifts  and  platform  misalignments,  and  also  those  errors  caused 
by  environmental  factors  such  as  ionospheric  delay  effects  and  gravi- 
tational deflections  and  anomaly.  The  vector  of  errors,  denoted  dx, 
is  modeled  as  the  difference  between  the  indicated  vector  of  state 
variables,  denoted  by  and  the  true  vector  of  state  variables,  denoted 


The  error  state  variables  are  modeled  by  first  order  linear  differ- 
ential equations  of  the  form 


TO 


-‘.y  ■L-'-  ■ 
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dx(t)  = F(t)^(t)  + G(t)w(t) 


where 


^(t)  is  the  vector  of  error  state  variables. 

F(t)  is  the  system  dynamics  matrix. 

G(t)  is  a gain  matrix. 

w(t)  is  a vector  of  white  noise  inputs. 

Since  the  error  state  variables  and  white  noises  are  assumed  Gaussian, 
they  can  be  completely  described  by; 

1)  the  mean  of  the  state  variable  vector; 


dx(t)  = £|^dx(t)j 


2)  the  covariance  of  the  state  variables; 


P(t)  = E^(^(t)  - ^(t))(^(t)  - dx(t) 


3)  the  mean  of  the  white  noises,  which  is  zero; 


E[w(t)]  = 0 


4)  the  covariance  kernel  of  the  white  noise  inputs; 


S(t^)^(t^  - tg)  = E[w(t^)w'^(t2)] 


where  5(t^  - t2)  is  the  delta  function. 


L _V  . 


Linear  measurements  are  available  of  certain  of  the  state  variables 


at  discrete  instants  in  time  and  are  modeled  to  be  of  the  form 


= H(tj^)^(tj^)  + v(tj^) 


(13 


generated  as  a difference  betv^een  INS  and  GPS  outputs,  where 
Z(tj^)  is  the  vector  of  difference  measurements, 

H(tj^)  is  the  measurement  matrix. 

v(tj^)  is  a vector  of  corruptive  v/hite  noise. 

The  discrete  time  corruptive  white  noise  is,  as  was  the  continuous 
time  w(t),  modeled  as  Gaussian  vrith  zero  mean  and  covariance  kernel 
described  as 


0 * ^ ^2 


(14 


It  is  also  assumed  that  the  system  white  noises  and  the  corruptive 
measurement  noises  are  uncorrelated  for  all  time. 


^l’*2 


(15 


For  notational  convenience  the  reference  to  the  independent  vari- 
able (t)  will  be  dropped  and  time  dependence  will  be  assumed  unless 
specifically  noted  otherwise. 

The  covariance  of  the  errors  must  satisfy  the  matrix  differential 
equation 


P = ft  + + gqg'^ 


(16 


I 


j.  between  sample  times  and  it  is  the  solution  to  this  equation  that  prop- 

agates  the  error  covariance  forward  in  time.  In  the  actual  implementa- 
tion G is  an  n X n identity  matrix,  and  thus,  in  the  above  differential 
T 

equation  GQG  will  be  denoted  by  £ only. 

The  Kalman  Filter  is  used  to  provide  an  estimate  of  the  errors 
being  committed  by  the  system.  This  estimate  is  corrected  at  discrete 
times  by  a Kalman  "gain"  weighted  residual,  - Hdx)  added  to  the 

estimate.  At  the  time  of  the  "update"  of  the  estimate,  a correction 
is  also  applied  to  the  system,  based  upon  the  filter  computed  gain. 

In  other  viords,  once  the  errors  are  estimated  and  measurements  are  taken 
to  correct  the  estimate,  the  Il’iU  is  driven  so  as  to  remove  the  esti- 
mated errors  from  it  also. 

Thus  the  propagation  of  the  estimated  error  covariance  forward 
10'  in  time  becomes  more  than  just  a solution  to  the  differential  equa- 

tion above,  but  must  also  include  the  measurement  update  equations, 
which  are  provided  by  Kalman  Filter  theory.  Propagation  eqmtions 
for  a covariance  analysis  are  presented  here. 

From  an  initial  condition  of 


the  covariance  is  propagated  forvfard  in  time,  and  also  between  measure- 
ment updates  by 


(19) 
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where  the  superscript  minuses  indicate  time  at  the  instant  of  taking  a 
measurement  and  just  prior  to  an  update.  The  covariance  is  updated 
then  by 


P - K^HP 


(20) 


where  the  plus  indicates  as  the  minus  signs  above  but  just  after  update 
incorporation.  The  General  Covariance  Analysis  Program  used  assumes 
complete  impulsive  control  over  those  states  which  are  updated.  This 
will  be  discussed  further  in  Appendix  B. 

Measurement  Underv/eighting 

While  the  above  equations  accurately  describe  the  basic  Kalman 
Filter  covariance  propagation  and  update  relationships,  one  change  has 
been  made,  in  accord  with  the  previous  studies  by  The  Analytic  Sciences 
Corporation  (TASC).  That  change  is  the  incorporation  of  measurement 
underweighting  (Ref  10«2.18-2.19) . This  is  a modification  of  the 
Kalman  gain  calcula.tion  as 


(21) 


where  u^  is  the  measurement  underweighting  factor  and  is  a scalar 
larger  than  one.  The  factor  used  is  equal  to  1.2  until  the  root  sum 
of  the  position  variances  (the  square  root  of  the  sum  of  the  position 
covariances,  P(l,l),  P(2,2)  and  P(3»3))»  is  less  than  1,000  feet. 
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When  this  occurs,  the  factor  is  set  to  one,  effectively  turning  off 
measurement  underweighting.  The  underweighting  remains  off  unless  the 
root  sum  goes  above  1,000  feet  once  again. 

Measurement  underweighting  is  used  to  prevent  the  filter  from  mak- 
ing sudden  changes  in  state  estimates  during  initial  transient  periods 
when  a new  sensor  is  acquired,  such  as  when  a satellite  set  changes. 
This  reduces  the  possibility  of  filter  divergence  and  will  also  reduce 
filter  sensitivity  to  a single  bad  measurement  (Ref  10»2.1S). 


» - • . • . - - • - ' • " • 


V,  BiHce  Shuttle  IMU  and  GFS  Truth  hodel 


Dsscri'oti  on 


The  design  Kalnan  Filter's  effect  cn  actual  systs"!!  pei'f orxance  and 
error  grouth  is  the  final  evaluation  of  its  effectiveness.  It  is  not 
possible,  or  even  desirable,  to  perform  this  check  on  the  actual  system 
during  the  many  iterations  on  the  'rarious  Kalman  Filter  parameters.  Be- 
cause of  this,  a simulation  model  of  the  actual  system  is  used,  which  is 
referred  to  as  the  trutli  model,  as  mentioned  in  Gnapter  IV,  the  truth 
model  consists  of  all  the  significanb  error  states  in  the  Space  Shuttle 
IKU  and  in  the  GPS  measurements. 

The  truth  model  used  in  this  study  is  a modification  of  the  truth 
models  used  in  earlier  studies  on  this  problem  (Ref  9*3«2-3»12,  10:4,9- 
4.I9,  and  11:2,13-2.14),  It  is  a modification  in  that  the  state  order 
is  somewhat  different,  certain  GPS  states  are  modeled  in  a slightly 
different  manner,  and  the  unique  portions  of  tlie  trajectory  used  plus 
the  ultimate  aim  of  the  study,  allowed  deletion  of  certain  states  used  in 
those  previous  studios.  The  deleted  states  are  for  external  aids  such  as 
TAGAN  and  taro  alulmeter,  as  viell  as  "pseudo  drag"  update  (Ref  95  4.4}, 


v.'hich  were  not  used  in  this  study. 

First,  the  error  states  are  presented,  followed  by  a discussion  of 
the  truth  model  system  dynamics  ratrix  and  the  associated  process  tioiss 
covariance  kernel.  I.'ext,  the  truth  model  data  base,  which  includes 
initial  conditions  and  modeling  '/alues  for  the  covariance  propagation 
equations,  vrlll  be  disciissed.  The  chapter  will  be  concluded  with  a 
presentetion  of  the  measurement  matrix  for  the  truth  mcdel  plus  the 
associated  corruptive  measurement  noise  covariance. 
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g^rror  States 

Table  I defines  the  error  states  and  number  of  error  sources  for  this 
study's  truth  model.  The  73-^^tate  truth  model  consists  of  36  Il-lb-rclated 
states,  3 states  for  fi-ravity  deflections  and  anomaly,  6 user  clock  states 
for  the  Crbitor  Master  Timing  Unit  (UTU),  and  28  satellite  error  states. 
These  last  23  states  are  actually  a quadruple  repetition  of  a 7-state 
model  for  each  of  the  four  satellites  used  during  each  update. 


Table  I:  Deorbit  Truth  Model  bi'ror  Btates 


Sect: or  States 


r or  Dou.rce 


■'n.  of  States 


1-3  Position  brrors 
4-6  Velocity  brrors 


7,3  User  Clock  Phase,  P’req. 
9-11  Platform  Misalignments 
12-14  ‘.ccelerometer  Biases 
I7-I7  Gvro  Bias  Drifts 


13-21  Itther  n':U  Clock  Mtatet 


22-24  .*Gcel.  Scale  Factors 

25-30  Accel.  Miisa ligaments 

(■•■issym.  and  Kon-orthog. ) 

31-33  Accel,  Non-Linearities 

34-36  jCravity  Deflections  and 
Anomaly 

37-42  Cyro  Mass  Unbalances 
(g-sensitive  drifts) 

43-45  Gyro  Aniscelasticlty 
(g^-censitive  drifts) 

46-73  Satellite  Lrror  States 


The  sections  referenced  in  Table  I refer  to  the  section  of  the  sys- 
tem dynamics  matrix  vrhich  helps  to  define  the  models  for  the  states  in 
that  section  (see  P’igure  3)* 


The  system  dynamics  matrix  v;ill  be  presented  in  three  sections,  as 
shonn  in  Figure  3*  and  listed  in  Table  I.  The  sections  break  the  error 
states  into  three  main  groupings.  First,  sections  Libeled  uhich 
relate  the  IKU-related  states,  including  the  gravity  model  errors,  vnll 
be  presented.  Next  will  be  a discussion  of  the  Crbiter  hTU  error  model 
deuicted  in  sections  3.,  followed  by  a descrintion  of  section  G,  which 
contains  the  GPS  satellite  error  models. 
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Figure  4.  Konzero  Blocks  Within  Partitions  A and  B of 
System  Dynamics  Matrix 

Figure  4 presents  an  expanded  vievi  of  the  partitions  of  Fig\ire  3 
la.beled  A , A , , A,,  B , B and  B,  . The  part  labeled  A , contains  the 

coupling  betvreen  error  states  in  the  A and  A,  sections,  ivhich  are  the 
IKU  related  states.  The  A portions  are  easiest  described  by  discussing 
the  nonzero  blocks  labeled  component  by  component.  Each  of  the 
blocks  will  be  presented  in  order  and  defined  by  their  reiational 
functions. 

Notation  developed  in  previous  studies  (Ref  9:]i',l-E,20)  vdll  be 
follov/ed  as  much  as  possible.  In  this  regard  the  following  definitions 
are  presented. 
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= Specific  Force  vector  in  P (platforiri)  franie 


is  the  3x3  transf oruation  matrix  to  go  from  the  P frame  to  the 
FGIF  (i)  (see  Chapter  II),  It  is  used  because  the  platform  misalignments, 
accelerometer  errors  and  gyro  drifts  are  defined  in  the  P frame,  while 
velocity  and  position  errors  are  defined  in  the  SGIF  frame. 

Block  S^  relates  velocity  errors  to  the  position  error  rates  and  is 
defined  as 


3^  = I -3  o 
-1  - 3x3 


Block  3^  is  defined  by 


LI  T 

RU  r RU 


v;here 


ji.  is  the  g!:ravitationa.l  constant 
I is  a 3x3  identity  matrix 

RU  is  the  3x1  position  vector  (in  EGIF,  see  Chapter  II) 

Block  3 relates  platform  misalignment  errors  to  the  velocity  error 

» 


rates  by 


vjhere 


3_  = T.  / F 
-3  -i/p  -p 
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(26) 


Block  5,.  relates  accelerometer  bias  errors  to  the  velocity  error  rates 

— S' 


ty 


S.  = I ^ ^ 

- 3x3 


(27) 


Block  Tv  coucles  accelerometer  scale  factor  errors  to  the  velocity 

“5 


error  rates  by 


S.  = T.  / F, 

-5  -i/p  -1 


(28) 


where 


0 

0 


0 

f. 


0 

0 


0 


(29) 


Block  ^ relates  accelerometer  misalifjnments  to  the  velocity  error 


rates  by 


^ = ^i/p  ^2 


(30) 
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where 
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£2-  0 


"3  0 


Block  relates  accelerometer  nonlinearitios  to  the  velocity  error 


rates  by 


where 


% ^ ^i/p  -3 


F = 

-3 


Block  ^ relates  gravity  deflections  and  anomaly  to  the  velocity  error 
rates  through  the  use  of  Xi/v»  developed  in  Chapter  II,  by 


So  = T.  / 
-8  -i/v 


Blocks  3,3.,  and  3.  vrill  be  described  in  the  next  section  of  this 
—a  — ab  — b 

chapter. 

Block  ^ relates  the  gyro  bias  errors  to  the  platforra  misalignment 
error  rates  by 


^ “ - 3x3 
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Block  couplec  the  jjyro  mass  unbalances  to  the  platform  misa3i‘ni' 
ment  error  rates  by 


-11  ^ 


vjnore 


0 


Block  3^^  relates  the  gyro  anisoelasticitles  to  the  platform  misalign- 


ment rates  by 


^11  = £5 


where 


Block  ^22’  located  along  the  diagonal,  relates  to  the  modeling  of 
the  gravity  deflections  and  ano:tialy  as  first  order  Gauss-iarkov  processes 
and  is  defined  as 


V V ^ • *.*• ' • * 


V,^  is  the  Crbiter’s  earth-relative  velocity  and 

du>  dv,  and  dw  are  altitude  dependent  correlation  distances  expressed 


in  a table  lookup  routine. 

The  process  noise  spectral  density  of  the  u'hite  noises  utilized  in 
the  modeljn^^  of  the  gravity  error  states  are  defined  by  the  desired 
steady  state  values  for  the  errors  modeled  here.  The  values  are  given  by 

a,  = (^1) 

(y.  is  the  standard  deviation  of  the  gravity  model  error.  The  effect  of 
3 

the  gravity  errors  decreases  v.dth  altitude  so  that  at  250,000  feet  the 
standard  deviations  have  been  attenuated  by  three  (Ref  9«3*'^-3»9»  70*21- 
3.22). 

Part  B is  defined  by  the  Crbiter  MTU  quartz  crystal  clock  error 
model,  shovm  in  Figure  5>  also  given  in  the  earlier  TA3G  studies  (Ref  10: 
2.23-2.25).  Figure  6 depicts  the  definition  of  the  3 blocks  and  Table  II 
defines  the  states. 
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Figure  3*  ^uertz  Crystal  Clock  Frror  Model 


In  Figure  5>  the  factors  labeled  and  are  constants  Khich  have 

-1  -U 

the  values  1.01  X 10  and  1 X 10  , respectively.  The  block  labeled 

"a"  is  the  shuttle  acceleration  along  the  major  thrust  axis  in  g's.  In 
this  study  "a"  was  approximated  by  the  acceleration  load  factor  along 
the  x-body  axis.  It  is  supplied  to  the  simulation  by  a table  looiiup  rou- 
tine as  a function  of  time  since  deorbit  ignition,  with  values  for  the 
table  from  a.  previous  study  (Ref  8j63,“2),  (Bee  Chapter  VII  for  discus- 
sion of  error  in  this  model.) 

The  process  noise  spectral  density  employed  in  the  simuJation  is 
presented  in  hcuation  42,  which  s’nov/s  the  linear  dependence  of  w^,  w^, 
and  upon  eacl'i  o', her,  i.e.,  their  correla Lio;). 
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0 
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1— • 

o 

0 

0 
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X 

10"^ 
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10"^ 

0 

21 

0 

0 

0 

0 

0 

0 

(42 


where 

p 

(ft/sec)  '/sec 

is  in  (ft/sec^)^/sec 
3 2 

is  in  (ft/sec"^)  /sec 

Part  G contains  the  satellite  states  with  each  satellite  modeled  by 
using  a modification  of  the  MTU  clock  model  plus  terms  for  range  bias 
errors  and  ionospheric  delay  effects.  This  vra.s  done  because  a model  for 
the  satellites  was  not  readily  available,  as  some  values  are  still  under 
test  and  evaluation  and  some  values  are  still  classified.  The  modifica- 
tions to  the  user  clock  model  were  made  so  that  the  numbers  used  would 
be  "close”  to  expected  performance  values  given  in  the  literature 
(Ref  1012.7-2. 8, 385. 2-5-12), 

The  modifications  led  to  the  use  of  a five-state  satellite  clock 
model.  This  is  the  same  configuration  as  for  the  MTU  model,  but  ’with 
the  acceleration  sensitivity  state  deleted,  with  two  additional  states 
added  to  this  to  account  for  range  bias  error  and  ionospheric  delay 
error.  The  dynamics  model  for  one  satellite  is  presented  in  Figure  7 
as  it  appears  in  the  system  dynamics  matrix.  The  clock  model  of 


Figure  5 minus  the  acceleration  sensitivity  state  is  an  accurate 
depiction  of  the  clock  portion  of  the  satellite  model.  Table  III 


defines  the 

states 

and 

the 

units 

of  the  states. 

Note 

that  Figure  7 could 

be  any  one 

of 

the  f 

’OUi' 

7x7 

blocks 

that  compose  block 

diagonal  G in 

Figure  3» 
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Glee:-:  Fiiase  L'rror 
Clock  Frequency  Error 
Aging  Biaa  Error 


iOCo 

ft/sec 

ft/sec*^ 


Flicker  Frequency  Error  ft/coc 

Random  Frequency  Rate  ft/sec"^ 

Error 


Ranre  Bias  Error 


Ionospheric  Delay  Error  I ft/sec 


Reiterating  that  although  the  satellite  models  do  not  rejjrcsent  a 
documented  model  for  the  proposed  cesium-standard  GPS  satellites,  the 
numhors  do  give  a fairly  representative  magnitude  of  errors  to  he 
expected  by  using  these  satellites.  The  magnitudes  for  the  modified 
clock  model's  initial  conditions  and  unite  noises  are  basically  scaled 


values  proportional  to  the  values  used  in  the  Orbiter  KTU  model.  The 
major  error  contributions  are  modeled  in  the  range  bias  error  states  and 
the  ionospheric  delay  states.  Equation  ^3  gives  the  process  noise  spectral 
density  used  in  each  one  ol  the  four  satellite  models. 


5.  E-09 


4.  E-O?  ij-.5  E-09 


0 0 


0 0 


0 0 


L ’^_0' 


4.p  E— 09 


4.  b-05  4.  E-07  0 0 

4.  E-07  4.5  E-09  0 0 


0 0 


0 0 


Note  that  the  scaling  factor  used  for  the  standard  deviations  from  trie 

_3 

KTU  clock  model  to  the  satellite  model  vra.s  10  , so  that  the  process 

noise's  spectral  density  differs  by  a factor  of  10  The  units  are  the 
same  as  for  (see  Equation  42). 


Truth  Nod  el  Data  "nase 

The  truth  model  data  base  includes  the  initial  conditions  on  the 
truth  model  states  as  initialized  covariances,  tho  associated  process 
noise  spectral  densities  where  applicable,  and  the  corruptive  measure- 
ment noise  strengths,  plus  any  time  constants  for  time  correlated  states. 
The  process  noise  spectral  densities  have  been  presented  ea.rlier. 

The  initial  conditions  have  been  derived  from  two  basic  sources, 
as  mentioned  in  the  introduction.  The  initial  conditions  are  presented 
in  Table  IV. 
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1 

Position-X""* 

5,417  feet 

2 

Position-Y^* 

2,844  feet 

3 

Fosition-Z^* 

9,332  feet 

4 

Velocity-X^* 

9.863  fps 

5 

Velocity-Y^* 

3.101  fps 

6 

Velocity-Z^* 

7.115  fps 

7 

ilTU  clock  phase 

100  feet 

O 

HTU  clock  freq. 

0.01  fps 

9-11 

Platform  misalign. 

I.377E-O3  rad/axis 

12-14 

Accel.  Biases 

I.6O7E-O3  fps^/axis 

15-17 

Gyro  Bias  Drifts 

I.697E-O6  rad/sec 

13 

I-iTU  Aging 

3.86S-O6  fps^ 

19 

Flicker  Freq.  Error 

20.2  fps 

20 

Random  Freq.  Err.  Rate 

1.37  fps 

21 

Clock  Accel.  Sensitiv. 

0.1  fpsVe 

22-24 

Accel.  Scale  Factors 

100  ppm/axis 

25-30 

Accel.  Misaligns. 

7.272E-O5  rad(s^  6 comp) 

31-33 

Accel.  i<’ on- Linearity 

I.OO9E-O7  fps^/g^ 

34 

Gravity  Anomaly 

function  of 

35 

Gravity  Deflection 

altitude,  see 

36 

Gravity  Deflection 

text,  Eq  4l 

37-42 

Gyro  Mass  Unbalance 

3.770E-09  rad/s/g 

43-45 

Gyro  Anisoelasticity 

3.77OE-O9  rad/s/g^ 

46-73 

Sat  clock  phase 

0.1  feet 

Sat  clock  freq. 

l.OE-04  fps 

j 

Sat  Aging 

3.36E-09  fps^ 

Flicker  Freq.  Error 

0.0201  fps 

Random  Freq.  Rate  Err. 

1.391:^-03  fps'^ 

Sat  Range  Bias 

3.5  feet 

Ionospheric  delay 

0.0024  fps 

*Kote  that  position  and  velocity  errors  have  strong  cross 
correlations  in  the  inertial  frame. 


v.V.V 


wiinsi’si 
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Notes  on  Batollite  'nodelinr  valnos 

The  method  used  in  choosin/f  the  values  for  the  clock  ncdel  portion 
of  the  satellite  modelin;v-  has  been  discussed  earlier.  The  filter  is 
driven  by  the  difference  betvieen  actual  measurements  and  estimated 
(predicted)  measurements.  The  satellite  contributed  errors  reflected 
in  this  difference  remain  to  be  discussed.  The  errors  are  modeled  in 
both  the  satellite  models  and  in  the  measurement  noise  vector,  V,  and 
thus  in  its  covariance  kernel,  R. 

As  explained  in  Chapter  III,  the  time  since  upload  of  the  satel- 
lites by  the  ground  station  at  Vandenberg  is  a factor  in  the  accuracy 
of  the  measurements  taken.  This  is  due  to  the  satellite  clock  drift 
and  satellite  ephemeris  error  (Ref  10:2,7-2.8).  Another  main  error 
contribution  is  the  ionospheric  delay.  The  ionospheric  delay  is  re- 
ported to  be  of  the  same  order  of  m3.gnitude  as  the  satellite  clock  er- 
rors (Hef  3s5“10). 

The  values  used  in  the  truth  model  assume  that  the  satellites 
have  not  been  recently  uploaded.  The  trajectory  of  the  Crbiter  during 
the  deorbit/reentry  portion  of  the  mission  simulated  occurs  over  almost 
one  fourth  of  the  earth’s  circumference  to  the  west  of  Vandenberg  AF3. 
It  is  thus  reasonable  to  expect  that  a large  portion  of  the  satellites 
used  would  be  those  that  had  not  yet  been  ui^loaded  (see  Chapter  III). 
This  may  not  be  as  true  during  the  post-blackout  tine  period,  but  this 
factor  was  not  taken  into  account  in  the  modeling. 

V/ith  this  in  mind,  the  values  for  the  range  bias  and  ionospheric 
delay  states  initial  conditions  were  chosen  to  reflect  values  of  3*5 
feet  a.nd  2,k  X 10  feet  per  second,  respectively. 
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Truth  Model  Measurement  Matrix  and  Measurement  Noise 


The  truth  model  measurement  matrix  reflects  the  three  areas  of  con- 
tribution to  error  in  the  difference  measurement,  which  are  IMU  errors, 
user  clock  errors,  and  satellite  errors.  An  advanced  design  GPS  receiv- 
er capable  of  simultaneous  measurement  processing  for  four  satellites 
has  been  assumed.  The  Kalman  Filter  incorporates  the  information  as 
eight  measurements;  four  each  of  position  and  four  each  of  velocity,  or 
one  of  each  from  each  of  the  four  satellites  in  use.  The  primary  meas- 
urement updating  interval  used  in  this  study  was  30*72  seconds  between 
measurement  updates.  This  value  was  chosen  arbitrarily,  buc  to  reflect 
an  interval  during  which  great  change  in  dynamics  would  not  be  encoun- 
tered, and  one  less  than  the  interval  in  which  the  satellite  constella- 
tion was  tested  for  best  geometry.  The  measurement  sequence  is 
assumed  not  to  start  until  just  prior  to  deorbit. 

As  noted  in  Chapter  IV,  the  measurements  (or  more  correctly,  the 
difference  measurements)  are  described  by 


Z = H dx  + V 


which  is  the  difference  between  actual  receiver  measurements  of  range 
and  range-rate,  and  the  predicted  measurement  using  INS  indicated 
values.  At  measurement  update  times  the  filter  corrects  its  estimate 
of  errors  by  adding  the  Kalman  gain  weighted  residual,  J^(Z  - H d^) , 
to  the  estimate. 

The  measurement  equation  is  composed  of  a set  of  four  pseudo-range 
and  four  pseudo-range-rate  equa,tions  which  are  derived  in  the  literature 
(Ref  2i  14-16,  73-^'^^,  and  7*4,50-4,53) • These  pseudo-range  equations 
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include  uuer  clock  phase  errors,  satellite  clock  phase  errors,  satellite 
ran-e  bias  errors,  and  satellite  epheneris  errors.  Terms  from  user 
clock  frequency  errors,  satellite  clock  frequency  errors,  user  and  satel- 
lite random  frequency  errors  and  an  ionospheric  delay  tern  are  included 
in  the  pseudo-range  rate  equation.  The  Crbiter's  position  and  velocity 
errors  are  also  included  in  these  equations  from  the  IMU  position  and 
velocity  states.  Satellite  ephemeris  errors  and  satellite  clock  errors 
are  actually  indistinguishable  froni  each  other  in  reception,  and  their 
effects  are  accounted  for  together.  Each  equation  also  has  an  associ- 
ated corruptive  vnnite  measurement  noise  torn,  part  of  V,  vfhose 
covariance  kernel  R is  included  in  the  Kalman  Filter  equations. 

The  measurement  matrix  is  of  the  form 


1 ' 

1 

|i 

1 

1 

1 

0 

1 0 1 

1 I 

1 

1 R 

1 

-r- 

1 

1 

1 

1 1 

' m 

2 

‘’1 

1 

1 1 

10.01 

1 

• 1 

’ • 

3 

ll 

1 

1 

•1oj- 

\ 

1 

1 T 

1 

1 

1 

0 

1 1 

1 

f D 

1 1 

! 

k- 

-s 

2 

1 

-r- 

1 * 

I 

— 

^.01| 

“ lr~\ 

1 

_L. 

1 

-1 

j_J-_ 

5 

,1 

I 

1 0 1 

1 

1 

1 

1 

rL. 

I 

0 

1 1 

1 

1 

1 D * 

► • 

6 

"’’'3 

1 

_L 

1 1 

10. oil 

----! — 

1 

- 4-- 

1 

I - 1 

— 1 - 



0 — — 

1 

V* 

< 

7 

ll 

1 

' 0 1 

1 

1 

1 1 

Hn 

1 

1 

0 

1 1 

1 

1 

1 1 D 

“3 

8 

1 

1 

|o.oi| 

1 

1 

1 

! 1 »: 

C'  > ; 
LN 


V/here  the  blocPis  K„  are  each  described  by 
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and  Ir.G.  in  the  unit  vector  comoonent  in  the  i direction  to  the  j 
J 

satellite  in  the  inertial  fra  Tie  (i). 

The  block  labeled  1 and  repeated  four  tines  in  is  given  by 


These  blocks  couple  the  appropriate  satellite  errors  to  either  positi 

or  velocity  states  during  meacurenent  updates. 

The  covariance  h of  the  corruptive  measurement  noise  is  diagonal 

with  values  used,  again,  from  earlier  studies  (Ref  lli2.14),  given  by 
2 

(12  ft)  for  each  of  the  four  range  error  equations  and  (0.05  ft/sec) 
for  the  four  range  rate  error  equations. 


i 

i 


] 


VI,  Kalman  Filter  Desit^n 


Basis  for  Kalman  Filter 

The  Kalman  Filter  is  based  upon  a reduced  order  model  of  the  system 
truth  model,  described  in  Chapter  V.  F.ajor  error  states  of  interest  are 
modeled  in  the  filter,  and  the  unmodeled  errors  from  the  system  are  taken 
into  account  by  changes  in  the  '..'Ey  the  major  errors  are  modeled.  This  is 
often  done  either  by  adding  vihite  noise  to  a state  model,  or  by  complete- 
ly changing  the  state  model,  such  as  by  changing  from  a bias  model  to  a 
first  order  Gauss-harkov  process.  The  filter  is  implemented  in  a reduced 
order,  or  suboptimal,  model  so  that  it  can  be  used  practically  for  online 
estimation  and  correction  of  the  error  states  of  interest. 

As  explained  in  Chapter  IV,  the  filter  makes  corrections  to  the  real 
world  I'.’S,  and  in  the  simulation  corrections  are  made  to  the  truth  model. 
Based  on  these  assumed  corrections,  the  Kalman  Filter  also  corrects  its 
ovm  estimate  of  the  error  states.  The  corrections  are  based  upon  meas- 
urements taken  from  an  aiding  source  external  to  the  IMU  at  discrete 
times,  which  in  this  study  is  the  GPf>  receiver.  The  source  should  have 
good  long  term  information  characteristics  to  offset  the  INC's  poor  long 
term  characteristics,  while  taking  ful3  advantage  of  the  INfi’s  good 
"short  term"  characteristics. 

In  the  INS,  the  dominant  mode  affecting  INS  characteristics  is  the 
84  minute  "Schuler  Pendulum"  oscillation  period.  Significant  errors  may 
arise  when  this  mode  is  observed,  especially  at  quarter  period  intervals. 
It  is  thus  over  a.  period  short  in  relation  to  these  intervals  that  the 
INS  exhibits  its  good  short  term  characteristics.  GPS,  on  the  other  hand, 
may  present  "noisy"  measurements  over  a short  term.  Taken  over  a period 
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long  enough  for  the  IKS  to  exhibit  error  growth  however,  the  character- 
istics are  quite  good. 

Filter  Description 

The  Kalman  Filter  error  states  used  in  this  study  are  a subset  of 
the  truth  model  states  described  in  Chapter  V.  Information  on  filter 
states  used  and  initial  design  values  is  based  largely  upon  earlier 
studies  (Ref  9i4.1-^.9,  10i4. 5-4.9,  and  11:2.10-2.12),  which  in  turn 
refer  to  NASA  filter  specification. 

Based  upon  the  suggested  designs  of  these  earlier  efforts,  an 
initial  design  of  a 17-state  filter  was  implemented.  Tuning  was  accom- 
plished on  this  17-sta,te  model.  Based  upon  insight  gained  from  the 
performance  of  this  design,  an  11-state  filter  and  an  C-state  filter 
were  also  implemented,  and  the  results  compared  to  the  performance  of 
the  17-state  design.  Since  this  study  examined  the  use  of  GPS  measure- 
ments both  prior  to  and  after  blackout  in  two  phases,  the  design  of 
the  filters  was  accomplished  in  the  region  prior  to  blackout,  and  then 
performance  was  examined  in  the  region  below  blackout.  Only  the  11- 
and  8-state  filters*  performances  were  examined  in  the  post-blackout 
region  because  of  the  11-state  filter's  performance  being  so  close  to 
the  17-state  filter's  performance  in  the  pre-blackout  time  frame. 

The  error  states  utilized  in  the  17-state  design  are  listed  in 
Table  V.  The  order  of  the  states  was  arranged  with  an  eye  toward 
removal  of  the  last  six  states  for  implementation  cf the  11-state  filter. 
Likewise,  removal  of  the  next  three  states  led  to  the  states  used  in 
the  8 state  design.  This  was  also  the  reason  that  states  7 and  8 were 
removed  from  states  18-21  in  the  truth  model  portion  dealing  with  the 


Table  V»  Kalman  Filter  Error  .States  and  Initial  Conditions 


State 

Description 

Standard  Deviation 

1 

Position-X^ 

1,330  feet 

2 

Position- 

6,339  feet 

3 

Position-Z^ 

20,507  feet 

4 

Velocity-X^ 

20,931  feet/sec 

5 

Velocity-Y^ 

6.525  feet/sec 

6 

Velocity-Z^ 

15*235  feet/sec 

7 

Clock  phase 

1,200  feet 

8 

Clock  Frequency 

0.05  feet/sec 

9-11 

Platform  Misalignments 

2.753E-03  rad 

12-14 

Accelerometer  Biases 

I.7OE-O3  feet/sec^ 

15-17 

Gyro  Bias  Drifts 

1,6972-07  rad /sec 

MTUf  The  particular  arrangements  of  which  states  to  drop  in  going  from 
one  design  to  the  next  were  again  based  upon  the  previously  mentioned 
TASG  studies  for  NASA  (Ref  9j4.1-4.4,  10:2.37-2.^0,  and  11:2.10-2.12). 
The  values  used  in  each  filter  design  are  discussed  in  the  sections 
dealing  with  each  individual  design.  Since  the  initial  conditions  did 
not  change  from  design  to  design  for  the  variances  of  the  error  states, 
they  are  presented  in  Table  V for  the  8-  and  11-state  designs  also. 

The  initial  conditions  for  the  position  and  velocity  states  were 
originally  given  in  a separate  coordinate  frame  not  used  in  this  study 
(Ref  11:2.11).  The  transformation  of  the  entire  6x6  initial  position- 
velocity  covariance,  including  cross  correlations,  was  performed  by 
computer  using  a transformation  matrix  provided  in  the  references 
(Ref  9*A,3).  The  resulting  complete  6x6  (lower  triangle  only)  initial 
covariance  for  this,  as  well  as  the  post-blackout  phase  is  presented 
in  Appendix  C. 

The  design  parameter  values  derived  from  the  efforts  of  this  study 
will  be  presented  next  along  with  an  explanation  of  why  they  were 


42 


chosen.  The  performances  of  each  of  the  three  filter  designs  will  be 
presented  and  compared  in  Chapter  VII,  As  mentioned  earlier,  the  designs 
used  are  based  largely  upon  studies  done  for  NA3A  by  TASG  (Ref 
and  Ili2,10-R.12),  The  suggested  filter  design  parameters  from  these 
studies  are  summarized  in  Appendix  G and  compared  there  to  the  final 
design  parameters  from  this  study. 
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17-F>tate  Filter 

The  design  of  the  17-state  filter  and  its  "tuning"  led  naturally 
to  the  designs  of  the  11-  and  8-state  filters.  General  comments  about 
the  method  of  design  and  tuning  of  the  17-state  filter,  especially  as 
it  relates  to  the  first  8 states,  will  apply  to  the  discussion  of  the 
other  two  designs  as  vjell. 

The  first  17  states  of  the  system  truth  model  are  modeled  in  the 
filter.  Where  the  models  used  for  these  17  filter  states  differ  from 
the  truth  model,  attempts  have  been  made  to  account  for  the  effects  of 
the  unmodeled  error  states  upon  these  17  states.  This  has  been  done 
in  three  general  v/ays. 

First,  the  initial  covariance  matrix  vras  changed  for  some  states 
to  reflect  a greater  uncertainty  in  the  knowledge  of  the  state's  condi- 
tion, This  is  important  if  the  filter  design  is  to  be  general,  as  the 
real  vjorld  error  states  may  not  be  as  small  as  the  values  that  have  been 
assumed  for  the  truth  model  in  this  study.  If  the  filter  begins  by 
underestimating  the  system  errors,  the  performance  may  not  be  such  that 
system  errors  are  bounded  and  divergence  of  system  errors  may  occur. 

Second,  the  effects  of  some  unmodeled  states  have  been  somewhat 
added  through  increases  in  the  strength  of  process  noise  on  some  states. 
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This  can  be  especially  effective  if  the  missing  state  can  be  approximated 
by  a bias. 

Finally,  in  two  of  the  groups  of  states,  the  accelerometer  biases 
and  the  platform  misalignments,  the  sta.te  models  have  been  changed  from 
zero  order  models  with  inputs  from  other  states,  to  first  order  Gauss- 
Karkov  process  models  (Bef  9:4.5), 

Specifically,  the  modeling  changes  for  the  various  states  are  as 
follows.  The  position  error  model  remains  the  same  as  the  truth  model 
except  for  the  increase  in  uncertainty  reflected  by  the  increase  in 
initial  conditions.  The  velocity  state's  error  model  is  changed  by  both 
an  increase  in  covariance  and  the  addition  of  process  noise.  The  proc- 
ess noise  helps  to  compensate  for  the  effects  of  truth  model  states  22 
through  36,  missing  in  the  filter.  Clock  state  S is  modified  so  that 
it  too  has  driving  corruptive  white  noise.  As  mentioned  above,  the 
model  for  the  pla.tform  misalignment  error  states  is  a first  order  Gauss- 
Markov  process  in  the  filter.  It  has  a very  large  time  constant  (not 
infinite)  with  driving  white  noise  inputs  from  the  gyro  bias  drift  rate 
errors.  The  accelerometer  bias  error  states  are  modeled  as  first  order 
Gauss-Markov  processes  too,  with  shorter  time  constants  and  driving 
white  noise  inputs.  The  models  for  the  gyro  bias  drift  error  states 
are  unchanged  from  the  truth  model. 

These  results  are  shown  in  Figure  8,  the  filter  dynamics  matrix, 
and  Table  VI,  which  gives  numerical  values  for  the  time  constants 
and  defines  the  values  for  the  process  noise  covariance  matrix, 
assumed  diagonal.  Note  in  Figure  8,  blocks  vjithout  entries  indicate 
zero-fill  for  tliat  block. 
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Note  that  state  6 has  a jifferent  value  of  noise  spectral  density 
than  states  4 and  5>  the  other  tvfo  velocity  states.  Since  the  velocity 
state  computations  are  performed  in  the  inertial  refei'ence  frame,  the 
geometry  of  the  Crbiter's  reentry  profile  can  give  some  insight  into  the 
reasons  for  having  a difference.  Note,  in  the  filter  dynamics  matrix 
(Fig  8),  and  also  in  the  system  dynamics  matrix  (Figs  3 and  4 and  Fqs 
25~33)»  that  the  velocity  errors  are  quite  dependent  upon  the  vehicle 
sensed  specific  force.  For  the  CFT-l  trajectory,  the  heading  is 
primarily  east-v/est.  Translated  to  the  inertial  frame,  it  can  he  said 
that  the  primary  velocity  components  of  the  Orbiter’s  inertial  velocity 
vector  v/ill  be  in  the  and  direction.  This  vfill  also  be  true  of 
the  Crbiter  acceleration  vector's  primary  components.  Thus,  larger 
errors  are  to  be  expected  on  the  and  velocity  channels.  Since 
smaller  errors  are  expected  on  the  velocity  channel,  to  estimate 
those  errors  better,  a smaller  ^ is  required  for  that  channel. 

Note  that  this  is  true  only  for  this  trajectory  and  would  certainly 
not  be  true  for  a polar  orbit  reentry.  However,  it  can  be  expected  that 
the  nature  of  the  Space  Shuttle  missions  will  be  quite  structured , and 
it  should  be  straightfon/a.rd  to  extend  this  idea  to  other  missions, 
V/hether  the  extra  accuracy  from  the  better  estimation  of  individual 
channel  velocity  errors  would  be  vrorth  the  effort  that  would  be  re- 
quired, is  not  within  the  scope  of  this  study. 

The  noise  spectral  density  values  for  states  7 and  8 are  quite  high 
compared  to  suggested  values  in  TASC  studies  (Ref  11j2,11)  (see  Appendix 
C).  The  '/alues  were  developed  iteratively  while  "tuning"  the  filter  to 
estimate  more  closely  the  observed  errors  in  corresponding  system  states. 
The  system  states  for  the  clock  were  affected  by  the  acceleration 


46 


sensitivity  state,  which  couples  in  the  component  of  lon*;ituclinal 
acceleration  on  the  shuttle  (see  truth  model  discussion,  Chapter  V). 

This  component  increases  gradually  prior  to  blackout,  until  at  blackout, 
its  effect  is  seen  in  a high  level  of  error  in  both  states  7 and  8.  The 
high  values  for  the  strength  of  the  process  noises  here  are  to  force  the 
filter  to  estimate  these  error  levels  more  correctly  over  a broad  range. 
This  is  discussed  more  fully  in  Chapter  VII  in  the  results  section. 

The  values  for  the  process  noise  spectral  density  for  states  9-11 
were  also  the  result  of  an  iterative  process.  They  are  related  to  the 
choice  of  initial  condition  for  the  misalignment  states  in  the  filter, 
chosen  to  be  twice  the  initial  misalignment  assumed  for  corresponding 
states  in  the  truth  model.  These  values,  along  with  the  presence  in 
the  filter  of  the  gyro  bias  states  allovjed  the  performance  of  the 
filter  to  closely  approximate  that  of  the  system. 

The  time  constant  value,  initial  conditions  and  noise  spectral 
density  for  states  12-14  are  suggested  values  from  an  earlier  TASC  study 
(Ref  in  which  the  modeling  of  these  states  was  addressed. 

Including  these  states,  as  well  as  states  15-17,  can  be  seen  laore  as  an 
aid  to  model  the  velocity  errors  and  platform  misalignments  correctly 
than  to  estimate  the  accclerometc'^  biases  and  gyro  bias  drifts  correctly. 
This  is  because  the  accelerometer  >iases,  along  with  the  pJatform  mis- 
alignments, represent  a major  error  contribution  to  the  velocity  states. 
The  gyro  bias  drifts  contribute  siruiif icantly  to  the  platform  misalign- 
ments and  thus,  airo  to  the  velocity  errors  (Ref  9»5*7)»  Thus,  by  add- 
ing the  last  six  states,  the  modeling  of  the  velocity  and  platform  mis- 
alignment states  have  been  rr;a.de  more  complete. 
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The  measurement  matrix  for  the  filter  is  identical  to  the  left  8x3 
portion  of  the  system  mefisurement  matrix.  It  does  not  change  from  filter 
design  to  filter  design.  Equation  48  describes  the  matrix. 


lif  = 


(48) 


V/here,  as  in  the  truth  model  measurement  matrix,  the  partitions  H „ . con- 
tain  vectors  to  each  of  the  4 satellites  being  used,  and  are  described 

by 


H„.  = 
-fj 


Ir  s . 
X J 


Ir 
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Ir  s , 

-g-J. 


Ir  s . 
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(49) 


and  Ir.s.  in  the  unit  vector  componet;t  in  the  i^*^  direction  to  the 
1 -i 

satellite,  given  in  the  inertial  frame. 

The  associated  corruptive  measurement  noise  matrix  does  however, 
change  from  filter  design  to  filter  design.  In  all  designs  it  is  8x3  and 
diagonal,  and  the  same  values  repeat  for  the  corruptive  noise  on  the 
filter  range  error  measurement  strengths  for  each  of  the  four  measure- 
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merits.  The  same  is  true  for  the  range  rate  measureraent  errors.  The  17- 
state  filter  measurement  noise  strength  values  are  given  in  the  following 
equations. 


Hfri  = 800  (ft)^  (50) 

Rfvi  = 0-05  (ft/sec)^  (51) 


11-State  Filter 

A 17-state  filter  presents  more  computational  burden  to  the  onboard 
computer  than  is  allowed  by  NASA  for  the  Kalman  Filter  function.  Thus, 
an  11-state  design  is,  at  present,  the  proposed  filter  to  bo  used  with 
GPS  as  an  external  aid  (Ref  11:2,10).  There  is  a significant  saving  in 
onboard  computer  workload  by  using  an  11-state  filter  instead  of  the  17- 
state  design.  This  savings  is  seen  in  the  numbers  of  storage  words,  and 
the  nuraber  of  adds,  multiplies,  divides,  etc.,  that  must  be  done  for  each 
iteration. 

The  11-state  design  parameters  were  developed  largely  from  insight 
gained  in  the  tuning  process  for  the  17-state  filter.  The  process  noise 
spectral  densities  for  this  design  are  listed  in  Table  VII.  The  initial 
conditions  and  time  constants  remained  the  same  as  those  given  for  the 
17-state  filter  (see  Tables  V and  Vl). 

Notice  that  the  process  noise  spectral  density  of  the  velocity  error 
states  remain  the  same  as  in  the  17-state  design.  There  does  not  seem  to 
be  a significant  need  for  any  change  in  this  parameter.  The  design 
philosophy  employed  is  conservative  in  estimating  the  errors  on  these 
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Table  VII.  ll-f>ta.te  Filter  Prooese  Noiee  3pectrai  Density 


ilta  te  Value 

1-3  0.0 

4,5  0.1076  (ft/sec)^/sec 

6 0.08  ( ft/sec )^/sec 

7 10.0  (ft)Vsec 

P See  Equation  52 

9-11  1*557  X 10  (rad)^/sec 


channels  and  also  on  the  clock  states.  This  is  due  to  observed  diver- 
gence tendencies  in  the  system  error  states.  The  level  of  noise 
strength  is  actually  higher  than  is  required  to  compensate  for  unmod- 
eled (in  the  filter)  error  states.  Thus,  the  level  of  noise  spectral 
density  used  in  both  the  11-  and  H-state  filters  for  the  velocity  chan- 
nels drives  the  velocity  errors  in  the  filter's  pre-update  estimates  to 
approximately  the  same  level  as  in  the  17-3tatc  design  (see  discussion  of 
results.  Chapter  VII ). 

The  noise  spectral  density  for  state  7 is  the  same  as  in  the  17- 
state  filter,  for  the  same  reasons.  It  vfas  desired  to  use  the  same  value 
as  in  the  17-state  design  for  state  8 in  the  3-  and  11-state  filter  clock 
models  also.  However,  this  va]ue,  in  conjunction  with  the  measurement 
noise  strengths  chosen,  resulted  in  numerical  difficulties  due  to  initial 
transient  conditions.  This  problem  was  overcome  by  changing  both  the 
process  noise  spectral  density  for  this  state  and  the  measurement  noise 
strengths  on  both  the  position  and  velocity  measurements  gradually  from 
values  which  did  not  cause  problems,  to  values  which  equaled  the  values 
used  in  the  17-state  filter.  These  changes  in  values,  as  functions  of 
time,  t,  since  the  initial  time,  t^,  are  reflected  in  the  following 
set  of  eo nations. 
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The  time  intervals  chosen  slowed  the  convergence  of  the  filter  and 


system  error  states  to  a steady  state  condition  sufficiently  so  that  the 


transient  problems  were  eliminatodi 


The  noise  spectral  density  for  the  misalignment  states  needed  slight 


increases  duo  to  the  removal  of  the  gyro  bias  error  states  which  directly 


contribute  to  the  corresponding  errors  in  the  system  model. 


In  the  post  blackout  time  period,  the  same  basic  design  is  still 


valid.  Because  error  buildup  during  the  time  of  radio  blackout  (and 


thus,  no  GPf>  measurements)  is  quite  severe  (Ref  11:2.13),  when  the  GPS 


measurement  sequence  is  reinitiated  there  is  again  the  possibility  of 


transient  problems.  Because  of  the  error  buildup,  the  initial  conditions 


are  somewhat  changed.  The  first  8 states  have  propagated  without  benefit 


of  the  high  accuracy  of  GP3  measurements  for  10  minutes.  Even  though 


drag  updating  is  assumed  to  be  used,  its  effect  is  only  to  slow  error 


growth  by  l/3  (Ref  9* 3.^-3. 3).  Finally,  the  longitudinal  component  of 
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Table  Vill,  Feet  Blackout  Initial  Conditions,  8-  and  11-State  Filters 


The  high  inirial  condition  on  state  8 reflects  error  buildup  prior 
to  blackout  as  v'ell  as  during  the  blackout  region.  Because  the  longitu- 
dinal acceleration  has  decreased  during  blackout  hovrever,  as  the  filter 
begins  estimating  the  error  on  this  state,  corrections  are  applied  to 
the  truth  model.  These  corrections  can  now  cause  a decrease  in  the  errors 
because  the  longitudinal  acceleration  effects,  coupled  through  the 
acceleration  sensitivity  state  of  the  truth  model  are  so  reduced.  In 
order  to  prevent  the  numerical  difficulties  during  the  transient 
conditions  mentioned  earlier,  and  to  prevent  the  filter  from  reducing 


its  estimate  of  the  error  too  rapidly  for  state  8 (v/hich  could  result  in 
underestimating  the  truth  model  error  state),  time  varying  values  for  the 
noise  spectral  density  for  state  3 and  the  measurement  noise  strengths 
were  again  applied.  These  values  are  summarized  in  the  following  set  of 
equations. 
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The  measurement  noise  strengths  used  allow  the  filter  to  heavily  vieight 
the  G?!5  measurements  and  to  apply  as  much  of  the  accuracy  available  in 
the  measurements  as  possible. 

S-Statc  Filter 

The  8-state  filter  needed  no  special  tuning  effort  to  achieve  per- 
formance comparable  to  that  of  the  11-state  design.  Limits  of  time  and 
resources  precluded  any  work  other  than  investigation  of  the  performance 
of  a filter  derived  using  the  design  parameters  of  the  11-state  filter. 


minus  the  three  misalignment  states.  Thus  the  discussion  of  the  11-state 
design  parameters  applies  for  this  section  as  well.  This  is  not  alto- 
gether unreasonable,  since  as  was  stated  earlier,  the  noise  spectral 
density  on  the  velocity  states  presents  a conservative  approach.  Be- 
cause of  this,  the  effect  of  the  misalignment  states  on  the  velocity 
channels  in  the  filter  is  minimized  in  the  11-state  design.  This  can  be 
seen  by  comparing  the  computer  plots  of  the  11-  and  8-state  designs  in 
Appendix  A,  for  both  regions  of  operation,  above  and  belovr  blackout. 


Baseline  Filter 

A final  run  las  accomplished  for  comparison  purposes  vdth  an  "opti- 
mal" baseline  filter.  This  vras  performed  by  operating  a filter  based 
upon  the  truth  model  within  the  Covariance  Analysis  structure.  Since 
this  is  a filter  design  which  completely  models  all  error  states,  it 
will  demonstrate  the  best  possible  performance  by  a Kalman  Filter  design. 

This  is  somewhat  misleading  because  the  GGAP  software  assumes 
complete  impulsive  control  of  the  states  being  estimated  by  the  filter, 
and  this  is,  of  course,  impossible.  However,  insight  can  be  gained 
into  the  performance  of  the  actual  suboptimal  design  filters  by  comparing 
their  performances  with  that  of  this  baseline  filter. 


VII,  Results  and  Conclusions 


Results 

This  section  presents  the  observed  filter  performance,  which  is  an 
examination  of  how  well  the  filter  was  able  to  control  truth  model  error 
growth  in  the  states  being  estimated.  The  computer-generated  plots  in 
Appendix  A will  be  the  focus  of  the  discussion,  since  they  allow  trends 
of  control  or  of  divergence  to  be  seen  at  a glance. 

As  briefly  noted  in  Chapter  V,  the  assumed  model  of  the  specific 
force  schedule  on  the  longitudinal  axis  of  the  Orbiter  for  the  clock 
model  was  modeled  incorrectly.  Due  to  a misinterpretation  of  a graph 
in  creating  the  table  lookup  routine  for  the  specific  force  schedule 
(Ref  8«63,82),  the  pre-blackout  region  had  a specific  force  schedule 
beginning  too  early  and  building  to  too  high  a value.  The  post-blackout 
region  did  not  have  enough  specific  force.  Unfortunately,  this  error 
was  not  discovered  until  the  study  was  complete.  The  results  of  this 
error  are  not  serious  enough  to  invalidate  the  basic  conclusions  of 
this  study,  but  trends  of  the  achieved  performance  values  are  in  error. 
The  basic  effect  has  been  to  make  the  pre-blackout  results  much  more 
pessimistic  than  they  should  be,  and  to  make  the  post-blackout  results 
much  too  optimistic.  The  performance  of  the  filter  is  not  in  question, 
but  the  performance  in  the  regions  for  which  the  results  are  presented 
needs  to  be  examined  in  the  light  of  the  above  information. 

An  attempt  is  being  made  to  correct  this  deficiency  in  the  study  by 
providing  a more  correct  schedule  of  acceleration  load  to  be  coupled  to 
the  clock  frequency  and  phase  error  states  by  the  acceleration  sensitiv- 
ity state.  The  results  of  this  effort  will  be  implemented  on  the  8- 
state  filter  design,  and  the  performance  of  this  design  will  be  examined 


under  these  more  correct  conditions.  When  complete,  the  results  will  be 
summarized  and  attached  to  this  report  as  an  addendum.  Meanwhile,  the 
results  of  this  study  will  be  presented  as  a function  of  the  assumed 
model.  See'  Appendix  D for  correction  addendum. 

The  results  are  presented  in  four  sections,  which  follows  the  same 
arrangement  as  the  plot  presentation  in  Appendix  A.  The  results  of  the 
"Baseline"  design,  or  optimal  73-2tate  filter  implementation,  will  be 
presented  first.  Next  the  performance  of  the  1? -state  filter  design  will 
be  examined,  including  a comparison  of  its  performances  under  two  updat- 
ing rates.  In  the  third  section,  the  performance  of  the  11-state  filter 
design  will  be  presented  for  both  the  pre-  and  post-blackout  regions  of 
employment.  Finally,  the  8-state  results  are  discussed  and  compared  to 
the  previous  designs  for  both  time  periods  of  operation.  Note  that  the 
17-Gtate  design  was  implemented  only  in  the  pre-blackout  region,  because 
the  11-state  filter  was  able  to  perform  so  well  in  comparison  in  that 
same  time  period. 

The  plots  present  a side  by  side  comparison  of  the  filter  indica- 
ted performance  of  the  system  error  states  being  estimated,  with  the 
actual  system  error  on  those  states.  For  each  filter  design,  the  per- 
formance of  states  one  through  eleven  were  plotted,  except  for  the  8- 
state  implementation,  for  which  only  the  first  eight  states  were  graphed. 

Optimal  Filter 

As  discussed  in  Appendix  3,  and  Chapter  VI,  an  implementation  of  all 
73-states  in  a filter  mechanization  within  the  software  structure  used  is 
unrealistic.  This  is  due  to  the  assumption  of  complete  impulsive  control 
at  update  intervals,  over  the  states  being  estimated.  However,  the 
implementation  does  give  some  insight  to  levels  of  performance  that  may 


56 


possibly  be  achieved  in  other  designs. 

The  most  important  state  in  the  truth  model  affecting  the  accuracy 
the  simulation  was  able  to  achieve  was  the  acceleration  sensitivity 
state  in  the  clock  model.  Its  effect  is  seen  largely  in  the  plots  of  the 
system  clock  frequency  error,  and  to  a lesser  extent  in  the  clock  phase 
error  plots. 

The  effect  of  complete  impulsive  control  in  the  optimal  filter  can 
clearly  be  seen  here  by  examining  the  plots  of  the  clock  phase  and  fre- 
quency error  (Figures  9 and  lO).  There  is  no  great  corruption  of  per- 
formance as  the  modeled  longitudinal  component  of  acceleration  increases. 
This  is  because  internal  states  in  the  clock  model  are  reset  at  update 
times.  This  is  really  not  physically  possible.  Thus,  to  the  extent  that 
this  occurs  here  and  in  the  satellite  model  portion  of  the  truth  model 
(which  states  are  obviously  impossible  to  reset),  this  portion  of  the 
simulation  is  not  realistic. 

The  results  here,  however,  can  be  used  to  identify  trends  and  to 
indicate  some  absolute  limits  in  accuracy.  For  example,  the  plots  of 
velocity  error  (Figures  14-16),  and  of  misalignment  angles  (Figures  17- 
19),  very  clearly  indicate  the  results  to  be  expected  as  real  world 
specific  force  increases  toward  the  end  of  the  pre-blackout  period. 

The  inability  to  estimate  the  misalignment  angles  because  of  lack  of 
sensed  specific  force  after  deorbit  thrust  cuts  off  is  also  evident. 

This  knowledge  is  useful  in  interpreting  the  results  of  the  actual  fil- 
ter design  implementations  in  the  follov/ing  sections. 


17 -State  Filter 

The  17-state  results  S'  rved  to  set  a level  of  performance  from  an 
initial  filter  design  that  v/as  expected  to  provide  the  "best"  results 
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of  the  designs  implemented.  Two  updating  periods  were  implemented;  one 
update  every  30.72  seconds  (determined  by  integration  step  size),  and 
one  update  every  15 *36  seconds. 

For  each  filter  design  (1?-,  11- , and  8-state),  the  performance  of 
the  user  clock  phase  and  frequency  error  states  was  due  to  virtually 
identical  model  implementations,  and  thus  the  discussion  in  this  section 
is  applicable  (in  the  pre-blackout  region)  to  other  sections  as  well. 

As  stated  in  the  previous  section,  the  most  important  factor  in  the 
accuracy  achievabl-;  was  the  performance  of  the  user  clock.  Specifically, 
the  response  of  the  user  clock  model  to  the  assumed  longitudinal  compo- 
nent of  acceleration  was  important  in  the  results  achieved  on  the  veloc- 
ity and  position  channels.  Recall  that  the  acceleration  sensitivity 
state  couples  this  component  of  acceleration  to  the  clock  frequency 
error  state,  which  in  turn  drives  the  clock  phase  error  state.  System 
performance  at  update  times  depends  on  these  two  states  and  thus  will 
vary  with  the  level  of  the  longitudinal  acceleration  component. 

Because  of  the  buildup  of  this  acceleration  component,  it  was  im- 
possible to  choose  time  invariant  filter  param.eters  that  would  allow  the 
filter  to  correctly  estimate  clock  frequency  over  the  entire  time  period. 
Thus  a compromise  estimation  level  was  chosen  (see  Figures  22-23) t Note 
that  the  filter  begins  by  "overestimating"  the  true  error.  This  is  to 
say  that  the  filter-indicated  error  is  larger  than  the  actual  error 
being  committed.  This  will  be  termed  a conservative  design  when  done 
intentionally.  As  the  acceleration  model  builds  up^  the  filter-indicated 
error  is  very  close  to  the  system  error  for  a time;  and  then  as  the  time 
period  approaches  blackout,  the  error  grows  and  the  filter  is  "under- 
estimating" the  actual  error,  or  indicating  less  error  than  is  actually 
being  committed.  Even  during  this  period  of  underestimation  the  filter 
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is  still  applying  positive  correction  to  the  system  error  at  update  peri- 
ods, and  is  bounding  the  error  quite  v/ell.  The  peak  error  observed  also 
corresponds  roughly  to  the  highest  level  of  acceleration  expected  during 
reentry  operations  on  the  longitudinal  axis  (Ref  8:63). 

It  is  evident  that  a filter  design  with  time  invariant  parameters 
for  the  clock  frequency  error  model  is  incapable  of  estimating  clock 
frequency  error  well  over  the  entire  range  of  acceleration  levels  expect- 
ed. The  level  chosen  here  is  a practical  compromise.  This  does,  how- 
ever, suggest  that  benefit  could  be  gained  with  an  acceleration  depend- 
ency for  oarameters,  notably  and  R-  for  the  velocity  states.  This 
will  be  discussed  further  in  the  Conclusions  section. 

The  clock  phase  error  plots  (Figures  20,  21)  show  that,  even  though 
a tendency  exists  for  the  errors  on  this  state  to  grow,  the  filter  is 
able  to  bound  the  error  growth  adequately.  The  filter's  estimation  is 
somewhat  conservative  to  allow  for  the  error  growth  tendency,  and  thus 
overestimates  the  system  error  over  most  of  the  range. 

The  effects  of  the  performance  of  these  two  states  is  seen  in  the 
filter  performance  on  the  position  and  velocity  error  states.  The  posi- 
tion states  performance  follows  very  closely  that  observed  on  the  clock 
phase  error  (see  Figures  24-29)  in  shape  and  in  magnitude.  Similarly, 
the  velocity  state's  performance  (see  Figures  30-35)  is  strongly  influ- 
enced by  the  user  clock  frequency  error  state.  An  example  of  how  the 
filter  copes  v;ith  this  influence  can  be  seen  in  the  plots  of  filter- 
indicated  and  true  error  for  the  y^  velocity  error  (Figures  32  and  33) » 
which  is  the  worst  case  performance  in  all  designs. 

The  estimation  is  initially  very  conservative.  The  system 'error 
states  on  the  velocity  channels  is  actually  corrupted  during  measure- 
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ment  updatee  at  certain  points  as  the  clock  frequency  error  grov.'s. 
filter*  of  course  does  not  have  knowledge  of  this  with  its  simplified 
model,  yet  continues  to  apply  a correction  to  the  velocity  errors  as  a 
"steady-state"  is  reached  at  each  level.  The  overall  effect  is  to 
bound  the  error.  The  rapid  growth  of  the  error  on  the  velocity  states  at 
the  end  of  the  period  is  due  to  the  increase  (actual)  in  specific  force 
as  the  Grbiter  enters  the  atmosphere.  The  filter  is  still  able  to 
apply  strong  corrections  at  update  periods  in  response  to  this  error 
growth. 

Another  effect  of  the  increase  in  specific  force  is  the  ability  of 
the  filter  to  estimate  misalignment  angles  more  effectively.  The  mis- 
alignment states  (see  Figures  36-^1)  performance  is  not  affected  much 
(if  at  all)  by  the  clock  states.  The  filter’s  estimate  of  errors  on 
these  states  is  very  conservative,  estimating  much  more  error  growth 
than  actually  occurs  in  the  system  states.  As  mentioned  in  the  discus- 
sion of  the  baseline  (optimal)  filter,  after  deorbit  thrust  is  cut  off 
4 minutes  after  deorbit  ignition)  until  specific  force  mounts  as 
atmospheric  drag  builds  the  filter  is  unable  to  estimate  misalignments 
to  a great  degree.  This  is  reflected  in  the  relatively  flat  area  on  the 
true  error  plots  for  these  states.  As  specific  force  grows,  however, 
the  filter  is  able  to  estimate  these  errors  more  effectively  and  is  thus 
able  to  apply  strong  correction  to  the  misalignments. 

Implementation  of  a faster  updating  rate  is  shown  in  the  next  series 
of  plots  (Figures  42-63),  3y  using  an  update  rate  vfhich  was  twice  as 
fast,  only  minimal  benefits  in  terms  of  accuracy  were  achieved,  as  seen 
by  comparing  these  plots  with  the  previous  set.  Table  IX  compares  the 
values  for  selected  system  error  states  at  the  next  to  final  update 
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before  blackout,  Pre-  and  post-update  values  are  presented  for  both  the 
15*36  second  and  the  30*72  second  update  periods.  V/hether  the  extra 
computational  load  would  be  worth  the  gains  in  accuracy  are  subjects  for 
further  study. 


Table  IX:  17-f>tate  Slow  and  Fast  Uodate  Rate  Comoarison 


System  (True)  Lrror  Std  Deviation 
17 -State  Filter,  t = t + 1351  sec 


Ufxjatc  = Ip 


2,  Position  ijrror  (feet)  | 23.7 
5,  Velocity  Error  (ft/sec )|  O.9I 


8,  Clock  Freq  Error  (ft/sec)  I 1,20 


Tj  / \ 

10,  Y*  Misalignment  (radians; 


8.16X10 


A final,  interesting  piece  of  information  that  can  be  gleaned  from 
these  performance  plots  is  the  effect  of  the  particular  satellite  con- 
stellation in  use  upon  accuracy.  This  effect  is  most  apparent  in  the 
plots  for  state  1,  X^  position  (Figures  24,  25).  Several  instances  of 
stronger  than  usual  corrections  to  system  error  are  immediately  apparent, 
especially  at  t 430  sec,  670  sec,  920  sec,  and  1350  sec.  These 
correspond  to  times  at  vfhich  the  satellite  selection  algorithm  switched 
to  a new  set  of  four  satellites,  based  on  GDOP  considerations.  This 
shows  dramatically  that  the  choice  of  satellites  can  have  a significant 
effect  upon  the  accuracy  that  can  be  achieved  using  GPS  measurements. 


11-State  Filter 

The  11-state  design  was  derived,  to  a large  degree,  from  the  design 
parameters  of  the  17-state  filter.  Much  of  its  performance  is  there- 


fore  very  similar  to  that  of  the  17-state  filter  (see  Figures  64-S5). 

The  major  insight  gained  is  that  the  effects  of  states  12-17  can  be 
adequately  compensated  by  state  noise  additions.  As  a matter  of  fact, 
the  plots  for  the  pre-blackout  region  show  that  better  estimations  of 
the  misalignment  errors  vras  achieved  in  this  design  (see  Figures  80-85)* 
This  resulted  in  the  11-state  design  providing  slightly  better  perform- 
ance on  the  velocity  channels  at  the  end  of  the  region,  when  the  mis- 
alignment states  vfere  rapidly  changing.  Other  differences  are  so 
quantitatively  small  that  they  will  not  be  summarized.  Because  of  this 
fact,  of  the  two  designs,  only  the  11-state  design  performance  was 
examined  in  the  post-blackout  region. 

In  the  post-blackout  region  (see  Figures  86-107)  the  modeling 
(incorrect,  see  chapter  opening  discussion)  of  the  longitudinal  compo- 
nent of  acceleration  on  the  Crbiter  as  significantly  decreasing,  allowed 
the  filter  to  apply  large  corrections  to  the  clock  frequency  and  phase 
errors  (see  Figures  80-83).  This  allovfed  substantial  improvement  over 
the  pre-blackout  region  performance.  Actually,  this  performance  is 
more  to  be  expected  during  the  pre-blackout  region  with  correct  modeling 
of  acceleration  inputs  to  the  clock.  In  contrast,  the  performance  of 
the  system  observed  during  the  simulation  of  last  half  of  the  pre- 
blackout region  would  be  the  expected  performance  in  the  post-blackout 
region. 

The  real  world  levels  of  specific  force,  sensed  by  the  D'U  allowed 
significant  corrections  to  be  applied  to  the  misalignment  states  (see 
Figures  102-107).  Their  effect  on  the  velocity  errors  was  reduced  to 
a great  degree,  Several  periods  at  vfhich  the  Shuttle  is  making  cor- 
rection turns  in  its  trajectory  become  apparent  in  their  effect  upon 


the  ability  of  the  filter  to  apply  la.rg'e  corrections  to  misalignment 
states.  This  is  due  to  increases  in  specific  force  in  all  directions 
during  these  turns. 

The  effect  of  the  reduction  in  the  misalignment  errors  may  not  be 
readily  apparent  until  the  11-state  results  are  compared  with  the  8- 
states  results.  This  v;ill  be  discussed  in  the  next  section. 


8-r>tate  Filter 

Surprisingly,  the  performiance  of  the  8-state  filter  is  very  close 
to  that  of  the  11-state  filter.  Before  discussion  of  the  plots  is  begun, 
a comparison  of  results  on  selected  states  between  the  11-state  and  8- 
state  designs  will  be  presented.  The  numbers  quoted  will  be  for  the 
error  state  values,  after  the  next  to  the  last  update  is  performed  for 
both  pre-blackout  and  post-blackout  regions.  These  are  given  in  Tables 
X and  XI. 
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Table  X:  Pre-Blackout  Comparison  of  11-  and  8-State  Filter  Performance 


System  (True)  Error  Std  Deviation  I 

t = t 

n 

+ 1351  sec 

11-Stat 

0 Filter 

Filter 

State 

+ ” 

+7 

2,  Position  Error  (feet) 

42.7 

24.8 

42.2 

23.1 

5>  Y^  Velocity  Error  (ft/sec) 

0.96 

1.13 

0.93 

1.13 

8,  Clock  Freq  Error  (ft/sec) 

1.94 

0.c<92 

1.94 

0.959 

10,  Y^  Misalignment  (radians) 

1.24X10"^ 

1.08X10”^ 

1.40X10"-' 

1.40X10"^ 
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Table  XI»  Poet-'^lackont  Go'''>arisor.  of  11-  and  P-State  Filter  Ferforrance 


System  (True)  ^rror  Std  De 

/iation 

t = t + 

{■ 

553  sec 

11-Sta  t 

a Filter 

' ate 

State 

X ' 

X 

. + 

2,  Position  Error  (feet) 

7.3 

6.9 

22.1 

6.9 

5»  Y“  Velocity  Error  (ft/sec) 

si 

1 

L 

0.046 

1.36 

0.04? 

8,  Clock  Freq  Error  (ft/sec) 

0.020 

_0^1Q 

0.020 

O^Olo  _ _ 

10,  Y^  Misalignment  (radians) 

1.06X10"^ 

1.01X10“’ 

l.33xio"-^ 

1.33X10"^ 

First  of  all,  notice  that  the  primary  differences  lie  in  the  magni- 
tude of  misalignment  angles  and  in  the  t”  values,  in  the  post-blackout 
region.  The  performances  after  updates,  at  t"*^,  are  virtually  identical. 
Also,  the  pre-blackout  region  results  indicated  very  little  difference 
existed  between  the  two  design  performances.  This  is  because  of  the 
real  lack  of  specific  force  inputs  to  the  IIW  and  thus  to  the  filter, 
and  thus  the  previously  discussed  inability  to  estimate  misalignments, 
until  the  end  of  the  period.  The  short  period  over  which  correction 
to  the  misalignments  does  occur  in  the  11-state  filter  is  not  long  enough 
for  the  effects  to  be  seen  in  the  velocity  and  position  errors.  In  a 
sense  then,  the  filter  implementations  are  almost  identical  for  much  of 
this  time  period. 

This  is  not  true  in  the  post-blackout  region.  A quick  comparison 
of  the  plots  for  the  S-state  filter  performance  (see  Figures  124-139) 
for  the  post-blackout  region  with  those  of  the  11-state  filter  (see 
Figures  36-107)  indicate  several  differences,  especially  in  the  velocity 
states.  The  plots  show,  and  Table  XI  confirms,  that  after  updates  have 
taken  place,  the  values  for  the  two  filters  are  very  close,  except  for 
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the  misalignment  states.  It  is  this  difference  in  misalignment  state 
error  values  that  is  responsible  for  the  large  difference  at  pre-update 
time  periods.  The  value  of  estimating  the  misalignment  angles  in  the 
post-blackout  region  can  be  seen  in  its  effect  of  reducing  the  rate  of 
error  growth  on  the  velocity  states,  and  thus  on  the  position  states 
also. 
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On  the  other  hand,  the  high  accuracy  of  the  GPS  measurements  pro- 
vides a post-update  level  of  accuracy  that  gives  the  8-state  filter  the 
ability  to  provide  an  average  level  of  error  on  the  position  and  veloc- 
ity states  which  is  quite  good.  With  the  corrections  to  the  clock  model 
acceleration  input  these  numbers  will  worsen  somewhat,  but  will  still 
be  acceptable. 

Conclusion 

There  are  two  main  conclusions  that  this  study  has  provided.  These 

a.rei 

• There  is  a need  to  model  filter  parameters  for  the  clock  frequency 
error  state  with  a sensitivity  to  longitudinal  axis  acceleration. 

• The  8-state  filter  design  can  perform  at  a level  which  v:ill  ade- 
quately bound  error  growth  in  the  IMU. 

These  two  conclusions  will  be  discussed  separately  in  the  next  few  para- 


gra  phs . 

Even  though  the  acceleration  modeling  for  the  input  to  the  clock 
sensitivity  state  is  in  error  regarding  the  time  at  which  it  is  applied, 
the  levels  modeled  will  be  encountered  at  other  points  in  the  mission 
profile.  Thus,  the  filter  will  be  called  upon  to  deal  with  essentially 
the  conditions  presented,  except  at  different  times  in  the  profile. 
Because  of  this,  the  filter  performance  on  the  clock  frequency  error 
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state  can  be  seen  to  be  a compi-onise.  If  the  filter  parameters, 

(the  filter  noise  spectral  density  for  state  number  8),  and  (the 
measurement  noise  strength  for  the  pseudo-range-rate  measurements)  vrere 
made  to  be  dependent  on  the  longitudinal  component  of  specific  force, 
the  filter  estimate  would  better  be  able  to  reflect  the  true  error  condi- 
tions. The  variation  around  the  estimation  level  chosen  for  the  design 
of  this  study  would  not  have  to  be  great  to  provide  much  improvement  in 
the  performance  observed. 

This  study  has  demonstrated  that  the  8-state  filter  is  able  to  per- 
form adequately  over  a wide  range  of  conditions.  The  effects  of  not 
estimating  the  platform  misalignments  during  the  short  time  period  after 
blackout  (and  also  during  blackout)  have  the  effect  of  slowing  conver- 
gence for  the  velocity  error  states,  and  through  them,  the  position 
error  states  also.  The  system  and  filter  are  not  as  sensitive  to  spe- 
cific force  inputs  as  when  the  misalignments  are  estimated.  These  fac- 
tors do  not  affect  performance  as  much  with  the  high  accuracy  GPS  meas- 
urements as  they  night  with  other,  less  accurate  external  aids. 

For  this  reason,  the  3-state  filter  should  be  given  serious 
consideration  as  the  design  choice  during  deorbit/reentry  operations, 
when  GPS  is  used  as  an  external  aid.  With  the  improvements  that  may 
be  possible  with  better  estimation  of  clock  frequency  errors,  the  8- 
state  design  should  provide  E\U  accuracy  well  within  required  limits. 
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This  appendix  will  present  the  computer  generated  plots  of  the  fil- 
ter performance.  There  will  be  seven  sections,  or  one  section  for  each 
filter  implementation.  First  will  be  the  baseline  or  73-state  optimal 
filter  plots.  Next  will  be  the  17-state  filter  with  30«72  second  update 
period,  followed  by  the  17-state  filter  with  a 15*36  second  update  peri- 
od. Fourth  will  be  the  11-state  filter  design  plots  for  the  pre-blackout 
time  period,  and  then  the  11-state  filter  results  in  the  post-blackout 
region.  Sixth  will  be  the  8-state  filter  performance  plots  in  the  pre- 
blackout region,  and  finally,  the  8-state  filter  design  plots  in  the 
post-blackout  region. 

Note  that  position  and  velocity  results  are  referenced  to  the  earth 
centered  inertial  frame  (i),  clock  errors  are  absolute  and  are  not 
referenced  to  any  frame,  and  misalignment  angles  are  referenced  to  the 
platform  (p)  frame.  For  each  section,  except  for  the  baseline  filter 
plots,  each  page  will  present  a side  by  side  comparison  of  filter  indi- 
cated error,  d$^,  for  a particular  state  with  the  true  system  error, 
dx^,  for  that  state.  The  baseline  filter  plots,  the  17-state  filter 
plots,  and  the  11-state  filter  plots  will  each  present  results  for  states 
one  through  eleven.  The  8-state  filter  results  will  be  given  by  plots 
for  states  one  through  eight. 

73-3tate  Optimal  Filter 

The  following  11  plots  (Figures  9-19)  show  the  performance  of  a 73- 


state  optimal  filter  implementation.  Only  one  plot  for  each  of  the 
first  11  states  is  shown  because  the  filter  and  system  would  have  the 


same  performance,  given  the  same  initial  conditions.  This  is  because  | 

the  filter  would  be  using  a near  perfect  model  of  the  system.  i 


tate  True  Error 
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17-State  Filter,  Slow  Ui^date  Rate 

The  following  section  presents  the  performance  plots  for  the  1?- 
State  Filter.  These  plots  are  for  states  one  through  eleven,  and 
represent  a simulation  with  a 30.72  second  measurement  update  period. 
The  first  time  point  for  which  data  vfas  taken  uses  120.55  seconds  as 
reference  or  t(o).  The  data  ends  at  assumed  blackout,  or  1482.71 
seconds. 
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ll-3tate  Filter,  ?ost-51ackout 

The  performance  of  the  ll-3tate  Filter  for  the  post-hlackout  portion 
of  the  simulation  is  plotted  in  this  section.  Initial  data  was  taken  at 
21,  seconds  and  final  data  was  taken  at  2704,  seconds.  The  measurement 
update  period  v.ras  30.72  seconds. 
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3-State  Filter,  Pre-Blackout 


The  8-State  Filter's  perfornance  during  the  pre-blackout  portion  of 
the  simulation  is  presented  by  the  following  plots.  Update  rate  and  time 
points  are  identical  to  those  of  the  11-State  Filter. 
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8-State  Filter,  Post-Blackout 

The  following  series  of  plots  shows  the  performance  of  the  3-State 
Filter  design  for  the  post-blackout  part  of  the  simulation.  Update 
rate  and  time  points  are  identical  to  those  of  the  11-State  Filter  in 
the  post-blackout  region. 
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Appendix  3 


Description  of  Computer  Programs  Used 

Software 

This  appendix  discusses  the  important  aspects  of  the  software  used 
in  this  study.  The  primary  program  was  The  General  Covariance  Analysis 
Program  (GCAP)  (Ref  6),  which  was  developed  by  personnel  of  the  Air 
Force  Avionics  laboratory  at  V/right-Patterson  Air  Force  Base,  Ohio  as 
a tool  to  aid  in  Kalman  Filter  design.  Because  of  the  special  nature 
of  the  simulation  of  the  GPS  satellite  constellations  required  in  the 
modeling,  the  software  was  modified  by  the  addition  of  several  sub- 
routines, The  added  software  was  taken  from  another  program,  developed 
by  personnel  from  a different  branch  of  the  Avionics  laboratory  to  aid 
in  GPS  receiver  design,  called  an  Integrated  GPS/lnertial  (IGI)  simula- 
tion, These  additional  subroutines  provided  the  simulation  of  the  24 
satellite  GPS  constellation,  the  selection  algorithm  which  examines 
visibility  and  "best"  satellite  configuration,  and  the  measurement 
algorithm  which  provides  Crbiter-to-satellite  vectors  for  each  of  the 
satellites  chosen. 

Each  of  these  two  programs  is  well  documented  (Ref  6,  15,  and  1?), 
and  thus,  neither  detailed  description  of  the  software  or  software 
listings  will  be  provided.  Instead,  discussion  will  be  limited  to 
general  descriptions  of  program  organization  special  features,  and  a 
discussion  of  the  satellite  visibility  and  selection  algorithm.  The 
appendix  will  be  closed  with  a short  discussion  of  the  implication  of 
using  a program  that  assumes  complete  impulsive  control  over  estimated 
states. 
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GCAP  was  designed  as  an  efficient  implementation  of  the  covariance 
analysis  equations.  It  takes  maximum  advantage  of  the  sparse  nature  of 
filter  and  system  dynamics  matrices,  measurement  matrices  and  noise 
covariance  matrices.  An  executive  routine  controls  the  operation  of 
the  program  from  initialization  to  final  output. 

Once  the  program  is  initialized,  a single  subroutine  for  integration, 
utilizing  a fourth  order  Runge-Kutta  integration  routine,  propagates  the 
covariances  of  the  filter  and  system  forward  in  time.  Elements  of  the 
dynamics  matrices,  noise  covariance  matrices,  and  measurement  matrices 
are  calculated  by  one  subroutine  for  the  filter  and  a separate  one  for 
the  system.  The  filter  and  system  each  have  separate  subroutine  groups 
which  perform  the  measurement  update  function,  called  from  the  integra- 
tion subroutine.  The  filter  subroutine  groups  controlling  updating 
compute  the  Kalman  gain  matrix,  1^,  and  correct  the  filter's  estimate 
of  the  error  states.  The  system  subroutine  groups,  on  the  other  hand, 
utilize  the  filter  computed  gain  matrix  to  perform  an  update  of  the  sys- 
tem error  states  which  v/ere  being  estimated  by  the  filter.  The  program 
is  completed  by  various  input  and  output  subroutines,  and  subroutines 
which  perform  the  mechanical  functions  of  matrix  multiplication,  in- 
version, transposition,  etc.,  specialized  for  both  sparse  and  full 
matrices.  Several  options  are  available  to  the  user,  dealing  with 
printirm  arj  plotting  functions,  and  the  frequency  with  which  each  is 
performed.  Th-'  program  performs  one  complete  cycle  of  integration,  up- 
date and  output  for  the  filter,  and  then  goes  back  and  does  the  same 
for  the  system,  after  which  the  entire  cycle  repeats.  Figure  pre- 
sents a generalized  functional  block  diagram  for  the  organization  of 
the  program. 


Figure  140.  Generalized  Functional  Block  Diagram  for  Modified  GGAP 


IGI 

The  subroutines  added  to  GGAP  from  the  Integrated  GPS/lnertial  sim- 
ulation software  perform  the  GPS  satellite  simulation  functions.  The 
primary  subroutine  of  this  group  generates  the  24  satellite  constella- 
tion and  propagates  it  forward  in  time.  This  subroutine  calls  another 
which  examines  the  total  set  of  satellites  to  decide  which  are  in  view 
for  the  Crbiter,  A satellite  is  in-view  if  it  is  10°  or  more  above 
the  Orbiter's  local  horizon.  After  the  "in-view"  set  of  satellites  is 


selected,  the  set  is  tested  iteratively  to  decide  which  four  satellites 
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from  the  set  will  provide  a geometrical  configuration  that  will  provide 
the  navigation  solution  with  the  best  accuracy. 

This  is  done  by  testing  for  the  best  "GDOP"  possible  from  the  in- 
view set.  GDOP  is  an  acronym  which  stands  for  "geometric  dilution  of 
precision".  It  is  defined  as  the  ratio  of  the  position  error  statistics 
to  the  range  error  standard  deviation  with  the  condition  that  all  sat- 
ellite range  measurement  errors  are  considered  uncorrelated  with  equal 
standard  deviations  (Ref  1:11).  GDCP  depends  only  on  the  geometry  of 
the  relative  positions  between  the  satellites  and  a user  and  does  not 
depend  on  the  range  error. 

Cnee  the  optimum  set  of  four  satellites  has  been  selected  it  is 
not  necessary  to  perform  this  test  again  until  the  relative  positions 
of  the  Orbiter  and  satellites  have  changed  by  some  amount  that  makes 
another  set  of  form  better.  The  program  provides  the  option  of  setting 
how  often  to  test  for  in-view  and  GDOP  conditions.  The  high  velocity 
of  the  Crbiter  during  initial  deorbit  operations  (25,000  fps)  required 
this  test  to  be  performed  fairly  often.  The  interval  was  set  in  this 
study  to  perform  this  test  every  61.4  seconds,  or  every  other  update 
period.  This  is  approximately  once  every  255  nautical  miles  of  Orbiter 
travel. 

The  measurement  subroutine  is  called  at  each  update  interval  during 
both  filter  and  system  update  operations.  In  turn,  it  calls  the  satel- 
lite generation/propagation  subroutine  to  an  entry  point  which  bypasses 
the  in- view  and  selection  portions  and  only  propagates  the  four  chosen 
satellites  forware  in  time.  After  this,  the  measurement  subroutine 
calculates  the  Orbiter  to  satellite  vectors  and  returns  them  to  the 
GCAP  updating  sequence.  The  set  of  subroutines  is  completed  by  a group 
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of  inversion,  multiplication,  and  transposition  subroutines. 

Impulsive  Control 

As  mentioned  earlier,  the  GGAP  implementation  assumes  impulsive 
control  over  the  estimated  states.  This  means  that  error  corrections, 
applied  only  to  software  in  the  filter  case,  are  assumed  to  be  applied 
instantaneously  to  a physical  condition,  in  the  case  of  the  system.  An 
example  would  be  position/velocity  indicators,  or  even  more  to  the 
point,  the  platform  misalignments.  These  physical  corrections  necessar- 
ily take  a finite  about  of  time.  Bias  levels  must  be  changed  and  phys- 
ical movement  must  take  place.  This  may  take  a significant  portion  of 
an  estimate  period  if  the  correction  is  large,  as  during  initial  tran- 
sients, By  the  time  the  correction  is  performed,  the  error  state  may 
have  propagated  forward  significantly  past  the  corrected  level  of  the 
error  state. 

Another  problem  is  the  assumption  by  GGAP  that  the  GPS  receiver 
clock  is  corrected  by  feedback.  This  requires  the  addition  of  an 
algorithm  which  will  send  estimated  errors  to  the  GPS  receiver  informa- 
tion processing  algorithms,  which  would  accomplish  this  feedback  func- 
tion. 

In  this  study,  the  changes  in  those  states  which  represent  phys- 
ical errors  are  small.  Thus,  total  impulsive  control  assumptions  are 
not  thought  to  present  significant  errors  in  the  analysis,  especially 
in  the  steady  state,  and  the  results  are  assumed  valid. 
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Appendix  C 

Filter  Parameter  Presentation  and  Comparison 

This  appendix  presents  values  for  filter  parameters  and  initial 
covariances  for  filter  and  truth  model  position  and  velocity  states.  No 
attempt  is  made  to  explain  differences  between  TA3G  suggested  values  and 
values  used  for  filter  parameters  in  this  study.  The  TASG  suggested 
values  are  from  a previous  study  for  NiiSA  (Ref  11j2,10-2,12) , 


Table  XIIj  Comparison  of  TaSC  Suggested  Filter  Process  Noise  and  Meas- 
urement Noise  Strengths  with  Final  Study  Values,  11-State 
Filter 
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Table  XIII » Full  Lower  Triangle  of  6x6  Initial  Covariance  Matrix,  Posi- 
tion and  Velocity  States,  at  Deorbit  I<^nition,  Filter 


6 


1.763X10'' 

-7.607x10'^  4.018X10'^ 
-2.653X10'° 


1.136X10'°  4.205x10'^^ 


2.462X10^  -1.098X10^  -3.338x10^  4.381X10^ 


-6.706X10' 


1.772x10^  -7.903x10^  -2.761x10 


2.991x10^  1.045x10^  -1.163x10^ 


4 2 

3.074x10 


4.258x10 


(72X10  2.321X10 


Table  XIV:  Full  Lovrer  Trianclo  of  6x6  Initial  Covariance  Matrix,  Posi- 
tion  and  Velocity  States,  at  Deorbit  I'>niticn,  Truth  Model 


1 

2 

3 
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2.934X10'' 

-1.^+6 1x10'^ 

3.039x10^ 

-4.965x10’^ 

2.506X10'^ 

8.708X10'^ 

-5.393x10^ 

2.745x10^ 

9.287x10^ 

9.728X10^ 

1.649x10^ 

-3.014X10^ 

-2.791x10^ 

-2.913x10^  9.615x10° 

-3.786X10^ 

1.918X10^ 

6.705x10^ 

6,896x10^  -2.079x10^  3.062x10^ 

Table  XV:  Full  Lower  Triangle  of  6x6  Initial  Covariance  Matrix,  Posi- 
tion and  Velocity  States,  at  Sxit-Blaokniit , Filter 


6 


1.17^I-X10  ' 


-3.072x10’^  2.575x10'^ 


8 


3.212x10''  -3.397x10  1.027x10 

1.213x10^  -5.412X10^  -1.891x10^ 


1.014X10' 


-3.304X10^  1.474X10^  5.150x10^  -2.136x10^ 


3.543x10^ 


3.731x10"^  -3.394x10'^  -1.361x10 


^ 3.108X10^  -5.957x10^ 


3.415x10^ 


Table  XVI:  Full  Lovfer  Triangle  of  6x6  Initial  Covariance  Matrix,  Posi- 
tion and  Velocity  States,  at  Mxit-31ackout,  Truth  Model 


6 


2.391x10' 

1.010X10^  6.626X10^ 

•2.994x10^  -1.136X10'^  7.04^«10^ 


5.599x10^  3.009x10^  -7.305x10^  2.078X10^ 


2.672x10^  1.674x10^  -2.346X10^*'  9.084X10^  5.246X10" 


-1.062x10^  -9.301X10  1.6‘^tXlO''  -2.761x10  -1.36‘^XIO  5.306x10" 
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Appendix  D 

Addendumt  Results  After  Model  Error  Correction 

Corrected  Specific  Force  Schedule 

As  mentioned  in  Chapter  VII,  an  error  was  made  during  the  original 
study  in  modeling  the  longitudinal  component  of  specific  force  on  the 
orbiter,  affecting  the  quartz  crystal  clock  frequency  error  of  the 
Master  Timing  Unit  (MTU)  model.  This  appendix  presents  the  results  of 
using  a correct  model  for  this  specific  force  schedule.  This  is  impor- 
tant, since  the  specific  force  model  provides  the  primary  effect  on  GPS 
clock  accuracy,  which  in  turn  influences  the  accuracy  that  can  be 
achieved  at  update  times  on  both  position  and  velocity  estimates.  The 
effects  are  depicted  by  plots  of  the  performance  of  the  8-state  filter 
designs  described  in  Chapter  VI  in  both  pre-  and  post-blackout  regions 
of  operation. 

The  corrected  specific  force  schedule  applied  to  the  clock  model  is 
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summarized  as  follows.  During  the  decrbit  ignition  period  there  is  an 
application  of  thrust  for  approximately  four  minutes,  followed  by  a 
period  of  free  fall  during  which  there  is  little  specific  force.  As  the 
trbiter  enters  the  atmosphere  specific  force  slowly  builds  in  response 
to  atmospheric  drag,  commencing  approximately  two  to  three  minutes  prior 
to  the  beginning  of  assumed  blackout.  At  exit  from  blackout,  the  longi- 
tudinal component  of  specific  force  has  increased  to  approximately  0.008 
g.  It  continues  to  build  until  reaching  a peak  of  approximately  0.56g, 
at  about  four  minutes  prior  to  landing,  and  then  rapidly  decreases.  At 
the  end  of  the  simulation,  when  the  Shuttle  is  at  20,000  feet  MSL,  the 
level  of  specific  force  is  0.3  g (Ref  Sj63). 

The  deorbit/reentry  portion  of  the  mission  is  one  of  high,  almost 
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constant  deceleration,  once  atmospheric  entry  is  accomplished.  This 
causes  a constant  increase  in  the  quartz  crystal  clock  model  frequency 
error,  and  also  in  the  clock  model  phase  error. 
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Filter  Kodiflcations 

Due  to  very  limited  time  and  resources  with  which  to  accomplish 
this  correction,  no  extensive  retiming  was  attempted.  At  the  same  time, 
however,  it  is  felt  that  only  minimal  gains  could  be  made  with  extensive 
retuning.  As  explained  in  Chapter  VII,  the  effect  of  the  error  was 
mainly  to  interchange  the  characteristics  of  the  pre-  and  post-blackout 
regions  of  operation.  Thus,  the  filters  were  only  modified  basically 
by  interchanging  values  of  noise  spectral  density  for  the  filter  clock 
frequency  error  state  and  for  the  filter  measurement  noise  on  the 
pseudo  range-rate  measurements  between  the  pre-  and  post-blackout 
regions  (see  ?g  5I  and  53  for  comparison). 

This  is  summarized  for  the  pre-blackout  region  as  follows 1 
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3.0  X 10 


3.0  X 10 


0.025 


6.0  X 10  ^ F (ft/sec)^/sec  < t - t < I50  sec  (58 


t - t <85  sec 
o 


t - t > 150  sec 
o 


0.009 


= < 0.035  > (ft/sec)  < t - t^  < 130  sec 


t - t <85  sec 

o — 


t - t > 130  sec 
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The  values  used  for  these  quantities  in  the  post-blackout  region  weret 
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In  all  cases  is  the  time  at  the  beginning  of  the  region  of  operation, 
either  pre-  or  post-blackout.  As  explained  in  Chapter  VI,  the  values 
•are  changed  gradually  to  prevent  filter  corrections  from  driving  truth 
model  covariances  negative  due  to  numerics  during  the  initial  transient 
period . 

It  is  net  intended  that  these  values  be  taken  as  final  values  for 
an  implementable  design.  They  should  rather  serve  to  indicate  possible 
levels  of  achievement  and  starting  points  for  final  refinement. 

Corrected  Filter  Performance 

The  filter  performance  will  be  discussed  in  reference  to  the  per- 
formance plots  for  the  pre-blackout  region  (Figures  141-156)  and  the 


post-blackout  region  (Figures  157-172).  Much  of  the  discussion  in 
Ciiapter  VII  is  applicable  here,  and  thus  only  major  differences  will  be 
pointed  out. 

In  the  pre-blackout  region,  the  largest  difference  is  that  the  sys- 
tem performance  is  muen  better  behaved  on  almost  all  states,  especially 
the  filter  clock  frequency  error  state  (see  Figures  111  and  144).  Cnly 
slight  improvement  is  noted  in  the  position  error  states.  This  suggests 
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a desirability  to  consider  using  a different  model  for  the  filter  posi- 
tion error  models  - perhaps  by  adding  pseudo  noise  to  these  states.  In 
the  present  form  the  achievable  accuracy  in  the  position  error  is  tied 
strictly  to  the  accuracy  achieved  in  the  velocity  error  states.  The 
velocity  error  states  show  much  improvement  over  the  earlier  results 
(see  Figures  121,  123,  15^»  and  I56) » due  mainly  to  the  improvement  in 
the  clock  frequency  error  state,  noted  earlier. 

In  the  post-blackout  region,  the  results  indicate  the  effect  of  the 
acceleration  sensitivity  in  the  quartz  crystal  clock  model.  The  overall 
interpretation  is  still  that  the  B-state  filter  performs  very  adequately 
in  bounding  error  growth  for  all  states,  supporting  the  conclusions 
drawn  earlier.  The  system  clock  frequency  error  especially  is  adequate- 
ly bounded,  so  that  its  effects  are  minimized  (see  Figure  I60).  A 
filter  utilizing  specific  force  dependent  parameters,  to  describe  the 
strengths  of  measurement  noise  and  process  noise  v;hich  affect  this  state 
should  provide  even  more  benefit,  however.  The  position  and  velocity 
error  states  show  only  slight  degradation  from  the  performances  of  the 
earlier  portions  of  the  study  (see  Figures  128-139  and  Figures  161-172). 
This  is  probably  due  to  the  influence  of  actual  specific  force  upon  the 
filter's  velocity  estimates.  This  enables  the  filter  to  estimate  the 
system  velocity  errors  better,  compared  to  the  pre-blackout  region  in 
which  specific  force  was  largely  absent. 

Conclusion 

This  appendix  supports  the  conclusions  drawh  in  Chapter  VII.  These 
are  that  the  3-state  filter  should  be  a viable  candidate  for  a GPS- 
aided  reentry  filter,  and  that  specific  force  dependent  filter  param- 
eters should  be  investigated  for  the  GPS  filter  clock  model. 


The  following  plots  (Figures  157-171)  depict  the  post-blackout 
performance  of  the  8-state  filter. 
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